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1.0  INTRODUCTION 


This  section  is  divided  into  two  parts.  In  the  first  part  we  discuss  the  motivation  for  this 
effort  and  in  the  second  part  we  briefly  discuss  fullerene  ion  thruster  work  being  carried  out  by 
others. 

1.1  Motivation  for  the  Present  Work 

The  motivation  for  the  present  work  stems  from  the  excellent  prospects  of  the  fullerene  ion 
thruster  (FIT)  to  become  the  preferred  electrostatic  propulsion  device  for  a  range  of  missions 
starting  from  low  thrust  attitude  control,  to  orbit  transfer  and  multimegawatt  interplanetary  flights. 
Simple  analysis  indicate  that  a  fullerene  Cgo  ion  thruster  could  provide  up  to  a  factor  of  30  increase 
in  thrust  density  over  the  current  atomic  fuels  of  choice  (Xe,  Hg)  as  well  as  a  simultaneous 
reduction  of  relative  losses  by  a  factor  of  5.5.  This  thrust  density  increase  would  result  in  a 
reduction  in  required  accelerator  grid  areas  to  a  size  that  is  economical  and  practical  to  construct 
even  for  megawatt  sized  thrusters. 

Space-charge  limitations  of  the  ion  acceleration  process  limit  the  thrust  density  of  Xe  or  Hg 
fueled  thrusters  to  the  extent  that  these  engines  are  presently  practical  for  limited  number  of 
missions.  Other  fuels  have  been  examined  to  increase  thrust  density  to  levels  that  would  lead  to  an 
attractive  thruster  with  wider  applicability.  Atomic  clusters  were  found  to  have  the  most  promise 
as  a  result  of  their  large  effective  molecular  weight.  Heretofore,  problems  with  cluster  stability 
during  acceleration  have  so  far  prevented  their  use  in  ion  thrusters.  However,  a  recently 
discovered  form  of  carbon,  termed  fullerenes,!1)  possess  a  wide  variety  of  unusual  properties. 
Their  high  molecular  weight,  extreme  stability  under  impact  plus  low  first  ionization  potential  make 
it  an  excellent  candidate  for  high  thrust  density  ion  thruster  fuel.  Below  we  briefly  substantiate  this 
assertion. 


As  mentioned  above,  an  intergrid  space  charge  limits  the  thrust  density  of  ion  thrusters. 
The  Child-Langmuir  space  charge  limited  current  equation  is 


j  = 


4^2 

- £ 

9  o 


/ 


.  m- 
\  1 


1/2 


V  3/2 
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where  Va  is  the  total  accelerating  potential,  d  is  the  grid  spacing  and  mj  is  the  ion  mass.  As  this 
equation  shows,  increasing  mi,  at  fixed  Va  tends  to  reduce  the  current  density.  However,  a  more 
likely  scenario  is  that  the  specific  impulse  (and  hence  the  ion  exit  velocity  (c),  would  be  kept 
constant,  since  an  optimum  value  of  c  exists  for  each  given  mission  and  power  supply.  In  that 
case,  the  accelerating  voltage  would  be  proportional  to  mj,  since 


V  =  — -i— - 
0  2e 

and  the  current  density  (Eq.  [1.1])  is  seen  to  increase  as  mj. 

The  grids  act  as  a  very  selective  filter,  extracting  ions  in  preference  to  neutrals  from  the 
partially  ionized  chamber  gas.  If  we  assume  that  only  the  ions  are  extracted,  the  mass  flow  rate  per 
unit  grid  area  is 
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Multiplying  times  c  gives  the  thrust  density  which  combined  with  Eq.  [1.1]  yields 


L-le(?L\2s 

A  9  °\  e  )  d2 


[1.2] 


This  shows  the  principal  advantages  of  using  high  mass  ions:  the  thrust  density  can  be 
made  higher  as  mj2,  thus  mitigating  one  the  principal  disadvantages  of  ion  thrusters. 

Extension  of  the  analysis  carried  out  above  lead  to  an  expression  relating  thruster  beam 
power  per  unit  grid  area  to  the  thrust  density  of  Eq.  [1.2].  This  expression  is 


EL.lt  L\c 

A  2  [AJ 


[1.3] 


As  an  example,  consider  an  interplanetary  mission  with  a  thruster  beam  power  of  10  MW, 
c  »  50,000  m/sec  and  grid  spacing  d  «  1  mm.  Equations  [1.2]  and  [1.3]  then  say  that  for  Xe,  the 
grid  area  is  18.3  m2  and  for  Cgo  it  reduces  to  0.61  m2.  The  later  is  practical,  the  former  is  not.  In 
fact,  it  may  be  advisable  to  reduce  the  thrust  density  for  other  reasons  and  still  maintain  practical 
size  grid.  Thus.  C^n  opens  new  doors  for  electrostatic  propulsion. 

Another  advantage  is  the  relative  reduction  of  power  losses.  These  losses  arise  from  a 
variety  of  sources  (ionization  work,  need  to  re-ionize  ions  that  strike  chamber  walls  before 
extraction,  electron  kinetic  energy  rendered  to  the  anode,  etc.).  In  general,  the  energy  loss  per  ion 
tends  to  be  independent  of  propellant  used,  since  it  is  closely  associated  with  the  ionization  energy, 
which  varies  relatively  little  (although  it  is  to  be  noted  that  C60  has  lower  ionization  energy  than  Hg 
and  Xe).  The  losses  per  ion  are  usually  quoted  as  a  voltage  AVl,  such  that  AVl  x  Ibeam  = 
Energy  Loss  Rate.  The  useful  work  ion  is  1/2  mjc2,  and  so  the  efficiency  is 


?1  = 
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[1.4] 


It  is  clear  that  higher  mj  leads  directly  to  higher  efficiency.  For  example,  using 
AVl  =  200  V  (a  conservative  estimate)  Fig.  1.1  shows  the  thruster  efficiencies  with  Xe,  Hg  and 
C60  used  as  propellants.  The  advantage  of  fullerene  is  most  striking  in  the  specific  impulse  range 
between  1000  and  2000  seconds,  where  many  space  missions  optimize,  and  where  the 
conventional  propellants  give  unacceptably  low  efficiencies  (as  well  as  low  thrust  densities,  (see 
Eq.  [1.2]).  The  dramatic  thrust  density  and  efficiency  improvements  with  C6o  fuel  are  shown  in 
Table  1.1. 
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Isp  (sec) 


Fig.  i.l  The  C60  Greatly  Improves  Ion  Thruster  Efficiency  in  the  Range  of  ISp  Where 
Most  Missions  Optimize  (1500  to  2500).  At  ISp  =  1500  the  C60  thruster  has 
double  the  efficiency  of  Xc  fueled  thruster. 
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TABLE  1.1 

Propellant  Comparison  for  c  =  50,000  m/sec,  d  =  1  mm,  AVl  =  200  volts 

VaXe/VaC^o  =  miXe/miC^o 


Propellant 

mi 

(g/mol) 

V 

(Volt) 

F/A 

(N/m2) 

Pjet/A 

(W/m2) 

0 

Xe 

131 

1700 

21.8 

5.46  x  105 

0.89 

C60 

720 

9340 

659 

1.65  x  107 

0.97 

Although  intergrid  voltage  approaching  10  kV  with  1  mm  spacing  may  not  be  achievable 
(in  test  chambers  at  p  ~  10"5  torr  the  peak  practical  voltage  may  be  2  kV  and  go  up  as  pressure  goes 
down),  Table  1.1  does  illustrate  the  Rillerene  ion  thruster  potential. 

1.2  Fullerene  Ion  Thruster  Work  by  Others 

The  potential  advantages  of  the  fullerene  ion  thruster  were  recognized  by  several  groups, 
some  of  which  carried  their  work  to  experimental  demonstration.  The  two  most  active  groups  are 
at  the  Jet  Propulsion  Laboratory(2"7)  and  Tokyo  Metropolitan  Institute  of  Technology^8-10) 

Numerical  modeling  was  performed  at  MITU1)  and  unpublished  exploratory  work  was 
performed  at  Hughes.  Our  own  initial  work  is  described  in  Refs.  12  and  13. 

At  JPL,  Leifer,  et  al(2-4>5)  focused  their  effort  on  investigation  of  the  fundamental  fullerene 
properties  with  fullerene  thermal  stability  being  the  principal  issue.  The  salient  results  were  that 
fullerenes  decompose  above  800°C  especially  in  contact  with  hot  metal  surfaces.  The  use  of  hot 
filament  cathodes  in  the  discharge  chamber  was  therefore  judged  undesirable  and  the  group 
constructed  laboratory  thrusters  with  RF  discharge.  The  experimental  work  was  primarily  carried 
out  by  Anderson,  et  al(3>6’7)  who  found  that  the  RF  discharge  established  on  Xe  quenches  when 
appreciable  fullerene  number  density  was  introduced  into  the  discharge  chamber.  This  was 
explained  by  known  electronegativity  of  the  fullerene  molecule  which  was  reported  to  scavenge 
low  energy  electrons  and  creating  negative  fullerene  ions.  It  was  ascertained  that  although  more 
power  could  be  delivered  to  the  plasma,  the  oscillatory  electron  motion  in  the  RF  discharge 
necessarily  creates  time  dependent  electron  energy  distribution  with  the  low  energy  electrons 
getting  captured  by  the  C60  and  therefore  applying  higher  RF  power  was  not  productive. 

Given  Leifer's  and  Anderson's,  et  al(2'6)  reasoning  for  selection  of  RF  discharge  over 
filament  cathode,  (i.e.,  avoid  hot  surfaces  inside  discharge  chamber)  the  next  step  would  be  to  use 
Electron  Cyclotron  Resonance  (ECR)  discharge.  This  however,  was  not  implemented  and  instead, 
Anderson,  et  al(7)  used  filament  cathode.  Using  this  system  he  found  that  up  to  2/3  of  the 
fullerenes  fragmented  and  remained  as  a  residue  inside  the  discharge  chamber.  The  unacceptably 
high  degree  of  fragmentation  was  attributed  to  electron  impact  ionization  process  which  is  in 
conflict  with  published  results  (Refs.  21-23  from  Ref.  7)  which  indicate  that  electron  impact 
fragmentation  occurs  at  energies  over  45  eV.  Since  the  discharge  voltage  was  about  40  to 
50  volts,  the  fragmentation  was  explained  by  long  residence  time  of  the  ions  in  the  discharge 
chamber  (~1  msec),  as  opposed  to  1  psec  in  the  experiments  where  the  45  eV  threshold  was 
measured.  Since  the  fullerene  ions  have  energy  dependent  half  life  for  dissociation,  Anderson,  et 
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alC7)  concluded  that  unless  techniques  for  much  shorter  residence  time  of  the  fullerene  ions  in  the 
thruster  discharge  chamber  are  found,  fullerenes  will  not  be  practical  propellant. 

The  Japanese  research)8"11)  lead  to  substantially  different  conclusions.  Their  work  also 
started  with  RF  discharged)  and  evolved  into  dc  discharge  using  filament  cathode.)9'11)  The  extent 
to  which  the  fullerenes  fragmented  and  how  much  if  any  of  the  sublimed  propellant  remained  as  a 
fragmented  residue  in  the  discharge  chamber  is  not  reported.  Their  best  results  are  36  mA  of  beam 
current,  extracted  with  1050  volts  intergrid  voltage,  at  12.2%  propellant  utilization  efficiency.  The 
beam  was  collected  on  a  target  plate  that  was  negatively  polarized  150  volts  below  the  accelerator 
grid.  The  discharge  chamber  diameter  was  approximately  11.5  cm  in  dia.  No  information  on  grid 
system  was  given.  Fragmentation  of  fullerenes  was  detected  in  the  beam  from  beam  spectral 
analysis  but  only  at  high  discharge  voltage  (60  volts  which  is  well  above  the  published  limit  of 
45  eV).  Reference  11  concludes  that  while  the  grid  system,  discharge  stability  and  plasma 
confinement  require  improvements,  the  fullerene  ion  thruster  is  viable. 
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2.0  IONIZATION  METHOD  &  THRUSTER  DESIGN  OPTIONS 


2.1  Background 

During  the  Phase  I  program  experiments  were  conducted  with  quartz  discharge  chamber 
with  2%  thoriated  tungsten  filament  cathode  at  its  center.C15)  The  filament  had  a  life  time  of  a  few 
minutes.  This  was  due  to  the  cathode  poisoning  caused  by  carbon  ion  bombardment  which  in 
some  cases  created  a  carbon  layer  as  much  as  0.05  mm  thick.  The  carbon  layer  then  became  the 
electron  emitter  which,  due  to  its  poor  emission  characteristic,  required  increased  temperature  of 
the  cathode  filament  eventually  causing  its  rapid  demise.  Thus,  cathode  selection  became  a 
significant  issue  with  broad  implications  to  the  design  of  the  overall  thruster  system. 

2.2  Discharge  Selection 

Fullerenes  fragment  to  graphite  when  exposed  to  temperatures  above  approximately  750°C. 
A  survey  of  literature  has  shown  no  elements  or  carbides  capable  of  thermionic  emission  at  or 
below  this  temperature. (16)  The  best  hollow  cathodes  tested  to  date  operate  at  about  1100°C. 
Thus,  operating  a  filament  cathode  or  a  hollow  cathode  on  C60  is  an  unlikely  proposition.  Without 
hollow  cathode  inside  the  discharge  chamber,  the  conventional  ionization  by  electron  bombardment 
must  be  replaced  by  some  electrodeless  discharge  such  as  RF  or  electron  cyclotron  resonance 
(ECR).  However,  the  RF  and  ECR  approach  requires  heavy,  inefficient  power  supplies  and 
discharge  chamber  walls  that  must  be  dielectric.  Since  some  C60  fragmentation  and  carbon 
deposition  on  the  discharge  cavity  walls  is  likely  to  occur,  it  could  lead  to  a  build-up  of  conductive 
layer,  inhibiting  the  RF  energy  transfer  and  shorting  the  discharge.  The  high  probability  of  this 
event  and  the  fact  that  the  external  neutralizer  is  likely  to  require  noble  gas  (at  least  for  long 
duration  missions),  lead  us  to  select  the  dc  discharge  method  with  conventional  hollow  cathodes 
operated  on  Ar  for  both  the  internal  discharge  and  the  external  neutralizer.  The  selection  process  is 
depicted  graphically  in  Fig.  2.1. 

Even  though  the  mass  flow  of  Ar  is  10%  or  less  of  the  Cgo  flow,  the  use  of  this  mixed 
propellant  with  very  large  differences  in  molecular  weight  may  have  some  undesirable 
consequences.  Among  them  is  system  complications  because  the  spacecraft  must  carry  dual 
propellant  and  the  possibility  of  larger  grid  erosion  due  to  the  C6o/Ar  interaction. 

The  increased  system  complexity  must  be  considered  relative  to  the  increased  power 
conditioning  mass,  and  inefficiency  if  the  dc  power  source  is  replaced  with  RF  or  ECR  power 
source.  The  preferred  approach  is  not  clear.  Initially  NASA  JPL  selected  the  RF  approach^2'5) 
which  was  a  factor  in  our  selection  of  the  dc  approach  because  it  offers  the  Air  Force  an 
opportunity  to  compare  the  two. 

The  C6o/Ar  interactions  were  considered  in  two  ways.  The  first  was  to  ask  what  could  be 
done  to  keep  the  two  species  apart  while  allowing  electrons  from  Ar  cathode  to  ionize  C60  vapor 
and  the  second  was  if  we  cannot  keep  them  apart  what  are  the  consequences,  primarily  to  the 
lifetime  of  the  grids?  These  issues  are  discussed  below. 

2.3  Alternative  DC  Discharge  Approaches 

A  conventional  ion  thruster  dc  discharge  chamber  has  a  coaxial  hollow  cathode  at  its 
upstream  end  that  operates  on  the  same  gas  as  the  main  discharge.  For  reasons  described  above, 
this  is  not  possible  for  the  fullerene  ion  thruster  (FIT).  We  therefore  considered  using  an  E-beam 
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and  also  placing  the  hollow  cathode  in  or  near  the  discharge  chamber  in  such  a  way  that  would 
allow  the  electrons  to  enter  the  C60  vapor  but  would  minimize  its  mixing  with  argon. 

Ionization  of  the  fullerene  vapor  with  E-beam  is  possible  but  a  disadvantage  is  the 
significant  complexity  of  the  E-beam  source,  high  electron  energy  and  low  current  density.  The 
current  density  can  be  increased  to  a  significant  level  only  by  increasing  the  electron  energy  near  or 
beyond  that  acceptable  for  fullerenes  which  is  about  50  eV.  Above  this  threshold  the  fullerenes  are 
likely  to  fragment  upon  C6o/electron  collision.  Because  of  this  difficulty  this  approach  was  at  least 
temporarily  abandoned. 

A  configuration  that  uses  hollow  cathode  as  the  source  of  electrons  for  C60  ionization  while 
limiting  the  mixing  of  the  hollow  cathode  gas  with  the  fullerenes  is  shown  in  Fig.  2.2.  It  is 
basically  an  Ar  (or  Xe)  ion  thruster  inside  a  C60  ion  thruster  that  share  a  common  hollow  cathode. 
This  approach  minimizes  C6o/Ar  mixing  within  the  discharge  chamber  and  offers  interesting 
operational  possibilities.  The  noble  gas  would  have  its  own  distinct  ion  optics  and  would  deliver 
high  ISp  and  low  thrust  while  C60  at  the  same  grid  voltage  would  deliver  higher  thrust  at  lower  Isp. 

This  may  be  most  attractive  for  large  ion  thrusters  with  power  levels  exceeding  10  kW  and 
grid  diameter  exceeding  50  cm  where  the  grid  structure  shown  in  Fig.  2.2  is  necessary  for 
mechanical  and  thermal  reasons  irrespective  of  the  nature  of  the  propellant. 

A  simpler  version  of  the  same  concept  is  to  eliminate  the  Ar  discharge  (anode  and  power 
supply)  and  the  Ar  grids  and  simply  exhaust  the  Ar  used  by  the  hollow  cathode. 

2.3.1  Negative  Ion  Beam 

Given  the  electronegative  nature  of  C^o,  it  is  natural  to  consider  construction  of  fullerene 
ion  thruster  where  the  Cgo  ions  are  negative  instead  of  the  conventional  positive  ions.  Negative 
C60  ions  can  be  created  by  flowing  the  fullerene  vapor  over  electronically  emissive  surface  which 
amounts  to  contact  ionization.  The  C60  acquires  low  energy  electron  (energy  approximately  equal 
to  the  emitter  work  function)  and  becomes  a  negative  ion.  The  negative  ions  can  then  be 
accelerated  just  like  the  positive  ions  except  for  reversed  grid  polarity. 

The  issue  then  becomes  how  to  neutralize  the  negative  beam  which  must  be  accomplished 
by  equal  beam  current  of  positive  charges.  One  possibility  would  be  to  simply  use  two  closely 
spaced  thrusters,  one  operating  on  negative  C6o  ions  and  the  other  a  conventional  noble  gas 
positive  ions.  Another  more  integrated  option  would  be  the  coaxial  dual  thruster  schematically 
shown  in  Fig.  2.2.  Electrons  from  the  central  hollow  cathode  would  impact  ionize  only  the  noble 
gas  and  openings  between  the  inner  chamber  and  outer  (C6o)  chamber  would  be  sealed.  The  C60 
would  be  charged  negatively  by  contact  ionization  as  described  above. 

The  fact  that  the  two  grid  systems  require  opposite  polarity  presents  an  engineering 
challenge.  The  good  feature  is  absence  of  external  cathode  neutralizer  which  is  not  needed  for  this 
dual  beam  thruster. 

2.4  Grid  Erosion  by  Mixed  Ar/C60  Propellant 

The  concept  shown  in  Fig.  2.2  was  evolved  because  of  our  concerns  over  unacceptable 
grid  erosion  caused  by  the  dual/mixed  propellant.  However,  before  embarking  down  this  path,  an 
assessment  of  the  magnitude  or  relative  grid  erosion  increase  was  needed.  Therefore,  ERC,  Inc. 
was  subcontracted  to  utilize  their  Plasma  Particle  Simulation  Code(17)  to  predict  what  is  likely  to 
happen  to  the  grids  due  to  direct  impingement  by  either  C60  or  Ar  and  due  to  charge  exchange 
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Fig.  2.2  Common  Hollow  Cathode  with  C60  and  Ar  Beams 
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collisions  with  corresponding  impingement  on  the  external  side  of  the  acceleration  grid.  The 


charge  exchange  collisions  that  were  considered  are: 

C6o+  (fast)  +  Ar  (slow)  -»  Cgo  (fast)  +  Ar+  (slow),  [2.1] 

Ar+  (fast)  +  Cgo  (slow)  -*  Ar  (fast)  +  C60+  (slow),  [2.2] 

C60+  (fast)  +  C60  (slow)  -*  C60  (fast)  +  C60+  (slow),  [2.3] 

Ar+  (fast)  +  Ar  (slow)  ->  Ar  (fast)  +  Ar+  (slow).  [2.4] 


However,  due  to  the  lack  of  collisional  cross-section  data,  only  the  fourth  charge  exchange 
reaction  was  actually  evaluated  and  resulted  in  negligible  molybdenum  grid  erosions  (8.5  x 
10"8  kg/hr  at  2000  volts  intergrid  voltage).  A  low  direct  impingement  of  C^o  was  also  predicted. 
The  calculations  were  performed  for  three  different  ion  optics  with  three  different  voltages  with  the 
Ar  mass  flow  being  10%  of  the  Cgo  mass  flow.  The  complete  ERC  report  is  reproduced  in 
Appendix  A. 

2.5  Penning  Ionization  Reaction 

Because  no  data  exists,  it  was  estimated  that  the  Penning  charge  transfer  reaction  of  excited 
Ar  and  Ar+  ions  at  thermal  energies  of  interest  has  a  rate  constant  greater  than  10‘9  cm^/molecule 
sec.(18)  This  is  a  very  fast  rate  which  means  that  at  2  microns  discharge  chamber  pressure,  which 
is  about  6  x  1013  Cgo/cm3  the  Ar+  reacts  to  make  C60+  in  about  17  x  10'6  sec  and  the  mean  free 
path  is  less  than  1  cm.  The  rate  constant  of  the  Penning  ionization  reaction  which  is  Ar  +  C60  -* 
Ar  +  C60+  has  not  been  measured  but  is  expected  to  be  greater  than  10‘10  cm^/mol  sec  in  which 
case  the  mean  free  path  at  2  microns  could  be  10  cm.  Thus  the  charged  Ar  ions  should  all  convert 
to  C60+  and  some  of  the  excited  Ar  should  produce  C60+.  However,  a  fraction  of  the  Ar+  may  be 
lost  on  the  walls  of  the  chamber. 

In  either  case,  the  concentration  of  Ar+  species  near  the  grid  should  be  small  and  therefore 
grid  erosion  due  to  direct  impingement  or  charge-exchange  collisions  should  be  small.  Most 
importantly,  the  ionized  and  accelerated  species  should  be  fullerenes  only. 

2.6  Section  Summary 

A  dc  discharge  with  hollow  cathode  operated  on  Ar  was  selected  for  fullerene  vapor 
ionization.  This  selection  was  based  on  (1)  assessment  of  electrodeless  discharge  techniques 
relative  to  the  dc  discharge;  (2)  survey  of  available  high  emissivity  elements  and  carbides; 
(3)  estimate  of  Ar+  concentration  near  the  grids;  and  (4)  estimate  of  dual  propellant  effects  on  grid 
erosion.  Although  this  approach  complicates  the  spacecraft  propellant  system  due  to  the  dual  fuel 
storage,  it  simplifies  the  neutralizer  design  and  operation  and  appears  to  have  greater  chance  of 
success  than  the  RF  discharge.  An  alternative  dual  beam  concept  was  presented  where  the  C60 
beam  uses  negative  ions  and  the  noble  gas  beam  uses  conventional  positive  ions  for  overall  beam 
neutrality.  Both  beams  of  course  produce  thrust  and  external  cathode  neutralizer  is  not  needed. 


10 


[  BUSEK  1 


BUSEK  Report  No.  036-3 


3.0  FULLERENE  STUDIES 


The  fullerene  vapor  pressure  data  that  were  obtained  during  the  Phase  I  investigation^5)  fell 
substantially  above  the  anticipated  values  (Fig.  3.1).  Because  vapor  pressure  (p)  and  temperature 
(T)  determines  the  propellant  mass  (m)  flow  which  is  a  crucial  variable  in  evaluating  space 
thrusters,  the  Phase  I  experiments  were  repeated. 

The  balanced  beam  apparatus  constructed  in  Phase  I  was  used  initially  to  measure  the  mass 
flow  real  time  and  which  was  then  related  to  the  vaporize  pressure  and  temperature  through  the 
equations  below. 


A Hs 


p  =  pre 

RT 

[3.1] 

m  =  pA^ 

[3.2] 

2tzRT 

where  A  is  the  orifice  flow  area,  Mw  is  molecular  weight,  R  is  the  universal  gas  constant,  AHg  is 
fullerene  sublimation  energy  and  Pr  is  a  reference  pressure.  AHs  and  Pr  are  given  in  Table  3.1 
below. 

Table  3.1  C60  and  C70  Vapor  Pressure  Related  Constants 


C60 

C70 

Pr  (mTorr) 

AHS  (cal/mole) 

Mw  (g/mole) 

4x  1010 
38,000 

720 

1.4  x  1012 
45,000 

840 

Prior  to  the  mass  flow  measurements,  all  fullerene  samples  are  first  "dried"  which  entails 
heating  the  batch  —  typically  few  grams  at  a  time  to  300°C  in  vacuum  (p  ~  1CH  torr)  for  about  24 
hours.  During  drying  the  sample  is  held  in  an  open  quartz  container  inside  a  vacuum  furnace. 

The  vapor  pressure  data  obtained  with  the  balanced  beam  apparatus  were  found  to  be 
inconsistent  with  poor  repeatability.  As  was  later  shown,  this  was  due  to  the  inconsistent  fullerene 
quality.  Furthermore,  the  fullerene  residue  —  (carbon  left  in  the  vaporizer  that  would  not  sublime) 
was  higher  than  previously  and  reached  as  much  as  30%  of  the  initial  load.  This  varied  from  batch 
to  batch  for  a  given  fullerene  supplier  and  especially  from  supplier  to  supplier.  High  residue  is  of 
course  unacceptable  (a  useless  mass  lifted  to  space)  so  it  was  decided  to  investigate  fullerene 
sublimation  and  purification  using  another  apparatus  that  is  simpler  and  removes  balanced  beam 
accuracy /thermal  drift  concerns. 

The  apparatus  is  shown  schematically  in  Fig.  3.2.  It  consists  of  two  separate  quartz  tubes, 
separated  by  a  replaceable  orifice.  The  tube  on  the  left  is  loaded  with  fullerenes  and  the  entire 
apparatus  instrumented  with  indicated  thermocouples  is  inserted  into  a  3"  dia.  tube  style  vacuum 
furnace.  The  section  to  the  left  of  the  orifice  is  located  in  the  hot  zone  of  the  furnace  and  serves  as 
the  vaporizer.  Most  of  the  section  to  the  right  of  the  choked  orifice  is  in  the  cold  zone  of  the 
furnace  and  serves  as  the  fullerene  condenser. 
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Typical  test  involves  loading  known  mass  of  fullerenes  into  the  vaporizer  section,  rapidly 
heating  to  specified  temperature  to  minimize  mass  loss  in  the  transient,  holding  temperature  for  an 
hour  or  longer  depending  on  the  test  objective,  followed  by  a  rapid  cool  down.  After  that  the 
fullerene  mass  in  the  condenser  and  the  mass  in  the  vaporizer  is  weighed.  The  vaporized  mass  and 
the  test  duration  are  then  used  to  calculate  m  and  finally  the  pressure  from  Eq.  [3.1]. 

The  results  of  these  tests  are  summarized  in  Table  3.2  and  plotted  in  Fig.  3.3. 

All  samples  purchased  from  three  different  suppliers  were  mixed  fullerenes,  typically  80% 
C60  and  20%  C70  The  suppliers  were  MER  Corp.,  Bucky  USA  and  Hoechst  AG  of  Germany. 
Of  these,  the  German  material  was  the  most  consistent  and  yielded  the  most  repeatable  data. 

The  quality  of  the  material  depends  on  the  production  method,  purification  method  and 
storage.  Of  these,  the  purification/extraction  solvents  are  the  most  important.  Typically,  solvents 
such  as  benzene,  toluene  or  CS2  are  used.  The  hydrocarbons  enter  the  C60/C70  cage  and  are  very 
difficult  to  remove.  Even  traces  can  be  responsible  for  polymerization  of  the  fullerenes  which  may 
remain  soluble  and  retain  the  basic  geodesic  structure  but  cannot  be  sublimed.  A  study  of  photo¬ 
illumination  of  C60  crystals  and  filmsU9)  lead  to  the  same  conclusion.  This  study  also  reports 
formation  of  polymerized  "skin"  on  the  C60  crystals  which  then  inhibits  free  sublimation  of  the 
interior  C60  thereby  lowering  vapor  pressure. 

The  conclusions  are: 

1)  Pre  qualification  of  fullerene  vendors  with  rigid  specifications  for  processing  and 
storage  is  be  required. 

2)  Drying  of  all  samples  is  required. 

3)  Purification  by  sublimation  is  probably  the  best  procedure  to  ensure  low  residue. 
(Most  thruster  tests  were  carried  out  only  with  in-house  sublimed  samples  using  the 
apparatus  in  Fig.  3.2.  For  large  quantities  the  vendor  supplying  the  fullerenes  should 
carry  out  sublimation  before  delivery). 

4)  Fullerenes  must  be  stored  under  inert  gas  in  dark  containers. 

5)  Direct  vapor  pressure  measurement  in  the  vaporizer  is  highly  desirable  but  no 
instrumentation  currently  exists  to  perform  this  measurement. 

6)  Sublimation  of  the  propellant  should  be  performed  at  as  low  a  temperature  as  possible 
and  on  as  small  a  fraction  of  the  total  propellant  as  practical.  This  preliminary 
conclusion  is  based  upon  uncertainties  of  Cgo  behavior  when  repeatedly  heated  and 
possible  photo-illumination  induced  "skin"  formation  by  hot  vaporizer  walls  (radiation 
from  metal  surfaces  at  or  above  600°C  is  partially  in  the  visible  spectrum). 
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4.0  LABORATORY  FULLERENE  ION  THRUSTER  (FIT)  DESIGN  AND 

ANALYSIS 


The  size  of  the  thruster  discharge  chamber  was  selected  to  be  approximately  15  cm  and  the 
active  grid  diameter  is  approximately  13  cm.  This  size  is  both  convenient  to  work  with  and  its 
performance  should  be  comparable  to  existing  fullerene  ion  thrusters  constructed  by  others!2-11) 
discussed  in  Section  1.2. 

This  section  contains  the  description  of  the  overall  design  and  performance  estimates  on 
which  the  design  is  based  followed  by  a  description  of  each  major  component  of  the  laboratory 
thruster  system. 

4.1  Overall  Design  and  Design  Analysis 

The  schematic  of  the  laboratory  fullerene  ion  thruster  designated  as  the  FIT- 13  is  shown  in 
Fig.  4.1.  Unlike  most  laboratory  thruster  experiments  the  propellant  storage  and  vaporization  unit 
is  inside  the  vacuum  tank  and  forms  structurally  integral  unit  with  the  thruster  itself.  We  are 
simultaneously  developing  both  the  thruster  and  the  propellant  feed  systems.  As  usual,  the  power 
supplies,  control  and  data  acquisition  system  are  located  outside  of  the  tank. 

The  sizing  and  design  analysis  of  the  FIT-13  system  with  two  different  grids  operating  at 
voltages  ranging  from  1000  to  3000  volts  are  shown  in  Tables  4.1  and  4.2.  The  equations  used  to 
make  these  predictions  are  summarized  in  Appendix  B. 

Table  4.1  shows  the  anticipated  thrust,  beam  current,  discharge  current  and  fullerene  mass 
flow  (along  with  other  parameters  of  interest)  for  a  set  of  grids  that  have  829  of  2  mm  dia.  holes, 
spaced  4  mm  apart.  This  set  has  been  fabricated  and  represent  the  initial  low  cost,  low 
performance  version. 

Table  4.2  shows  the  important  design  parameters  for  the  same  thruster  with  more  advanced 
grids.  The  thruster  system  has  been  designed  to  accommodate  this  new  higher  mass  flow,  beam 
and  discharge  current  intergrid  voltage,  etc.,  listed  in  this  table.  This  set  of  grids  may  be  fabricated 
in  future  programs. 

4.2  Propellant  Storage  and  Vaporization 

A  spacecraft  with  fullerene  ion  thrusters  performing  a  low  to  GEO  orbit  transfer  may  have 
on  board  hundreds  of  kilograms  of  fullerenes.  How  to  store  the  propellants  and  how  to  feed  it  into 
the  vaporizer  is  a  challenging  engineering  problem.  Several  approaches  can  be  considered  which 
defer  depending  on  the  spacecraft  mission  and  hence,  the  amount  stored  . 

Fundamentally  they  can  be  divided  into:  (1)  heating  the  entire  fuel  storage  in  which  case  the 
storage  tank  becomes  the  vaporizer;  or  (2)  heating  small  portions  of  the  propellant  in  a  small 
vaporizer  that  is  fed  incrementally  or  continuously  from  the  storage  tank. 

Both  approaches  have  drawbacks  and  advantages.  The  obvious  disadvantage  of  the  first 
approach  is  that  one  is  heating  large  volume  (of  the  order  of  0.6  m3  for  1000  kg  of  propellant) 
which  implies  large  heat  loss  and  large  thermal  inertia. 
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Table  4.1 

13cm  Laboratory  Fullerene  Thruster 


INPUT/CASE  NO 

1 

2 

3 

4 

!  5 

Dact.(m)  grid  active  dia 

0.13 

0.13 

0.13 

0.13 

0.13 

ds(m)  screen  grid  hole  dia. 

0.002 

0.002 

0.002 

0.002 

0.002 

da(m)  acc.  grid  hole  dia. 

0.002 

0.002 

0.002 

!  0.002 

i  0.002 

p(m)  grid  hole  pitch 

j  0.004 

0.004 

0.004 

0.004 

!  0.004 

lg(m)  grid  spacing 

0.0020 

0.0020 

0.0020 

!  0.0020 

j  0.0020 

Vt(Volt)  total  inter,  grid  voltage 

1000 

1500 

2000 

2500 

3000 

Vn(volts)  net  voltage=mi*cA2/e 

500 

1000 

1500 

2000 

2500 

eatu(-)  propellant  utilization 

0.8 

0.8 

0.8 

0.8 

0.8 

Tn  (k)  guess  at  dis.  cham.  Tvapr-ITERATE 

662 

676 

687 

695 

702 

epsbi(volts)  loss  per  beam  ion 

500 

500 

500 

500 

500 

Vdis.  (Volts)  quess  discharge  voltage 

i  30 

30 

30 

30 

30 

etabase(-),  base  eff. 

i  0.8 

0.8 

0.8 

0.8 

0.8 

Mw(kg/kgmol)  propellt  mol.  wt. 

720 

720 

720 

!  720 

720 

dstar  vap(mm)  vaporizer  orifice  dia. 

5.8 

5.8 

5.8 

5.8 

5.8 

Mwcathode(kg/kgmole)  c.  gas  mole  weight 

39.9 

39.9 

39.9 

39.9 

39.9 

fcn=(n  of  Ar)/(n  of  C60)  cathode  flow  fract. 

3.62 

2 

1.31 

0.94 

0.72 

(vary  fen  to  get  at  most  3sccm  cathode  flow) 

OUTPUT -THRUSTER 

fo(-)  open  area  fraction 

0.196 

0.196 

0.196 

0.196 

0.196 

Aat(m2)  surface  area  of  grid 

1.327E-02 

1.327E-02 

1.327E-02 

1.327E-02 

1.327E-02 

Aao(m2)  open  area  of  acc.  grid 

2.604E-03 

2.604E-03 

2.604E-03 

2.604E-03 

2.604E-03 

N  (-)numbr  of  holes 

829 

829 

829 

829 

829 

mi(kg)  mass  of  prop,  ion 

1.195E-24 

1.195E-24 

1.195E-24 

1.195E-24 

1 .1 95E-24 

ci  (m/sec)  ion  exit  speed 

1.158E+04 

1.637E+04 

2.005E+04 

2.31 5E+04 

2.589E+04 

le(m)  effective  grid  spacing 

2.236E-03 

2.236E-03 

2.236E-03 

2.236E-03 

2.236E-03 

J(A/m2)  beam  current  density 

1.289E+01 

2.367E+01 

3.645E+01 

5.093E+01 

6.695E+01 

Ibeam  (A)  beam  current 

0.034 

0.062 

0.095 

0.133 

0.174 

ni(1/m3)  ion  number  density 

6.948E+1 5 

9.025E+1 5 

1.135E+16 

1.373E+16 

1.614E+16 

cn(m/sec)  speed  of  neutrals 

124 

125 

126 

127 

127 

ni/nn(-) 

0.0427 

0.0305 

0.0251 

0.0219 

0.0197 

nn(1/m3)  neutral  number  density 

1.626E+17 

2.957E+1 7 

4.516E+17 

6.275E+1 7 

8.207E+1  7 

ntot(1/m3)  tot  number  density 

1.696E+17 

3.047E+1 7 

4.629E+1 7 

6.412E+17 

8.368E+1  7 

ntot(1/m3)  based  on  Tdis.  chamb.  input 
ntot(=R42)/ntot(=R43)iteterate  to  =1 

Tdis.  cavity  (degC)  based  ntot 

1.650E+17 

1.028E+00 

389 

2.940E+1 7 

1  036E+00 

403 

4.551E+17 

1 .01 7E  i-OO 

414 

6.199E+1  7 
1.034E+00 

422 

8.076E+1  7 

1.036E+00 

429 

p(N/m2)  in  dis.  cavity  based  on  ntot 

1.549E-03 

2.843E-03 

4.389E-03 

6.150E-03 

8.1 07E-03 

p(Torr)  in  dis.  cavity  based  on  ntot 

1.154E-05 

2.118E-05 

3.270E-05 

4.582E-05 

6.040E-05 

F(N)  thrust  from  eq.4 

0.0029 

0.0075 

0.0142 

0.0229 

0.0337 

F(grams))  thrust  from  eq.4 

0.30 

0.77 

1.45 

2.34 

3.44 

mdot(kg/sec)  mass  flow  eq.5 

3.1 29E-07 

5.749E-07 

8.851  E-07 

1.237E-06 

1.626E-06 

mdot(mg/sec)  mass  flow  eq.5 

0.313 

0.575 

0.885 

1.237 

1.626 

Isp(sec)  specific  impls 

947 

1338 

1638 

1891 

2114 

etat(-)  thrust  efficiency 

0.40 

0.53 

0.60 

0.64 

0.67 

Pttl=FA2/(2*mdot*eta) 

33.73 

92.80 

190.37 

332.43 

524.24 

Pdis=Pttl-Pbeam=Pttl-lbeam*Vnet 

16.95 

31.17 

48.04 

67.21 

88.45 

ldis.(amps)=Pdis/Vdis 

0.57 

1.04 

1.60 

2.24 

2.95 

OUTPUT -VAPORIZER 

guess  Tv(degC) 

625 

653 

673 

689 

703 

Pc60(mTorr) 

2.23E+01 

4.25E+01 

6.59E+01 

9.22E+01 

1.23E+02 

Pc60(N/m2) 

3.00E+00 

5.71  E+00 

8.84E+00 

1.24E+01 

1.65E+01 

mdot  c60  (mg/sec) 
mdot/mdot(from  row#52) 

(iterate  Tv  to  get  Mdot/mdot=1) 

3.10E-01 

0.99 

5.82E-01 

1.01 

8.91  E 

:-oi  j 
.01 

| 

1 

241 

E+00 

1.00 

1.63E+00 

1.00 

OUTPUT  -  cathode  max  flow 

Nc.gas(1/m3)  cathode  gas  number  densty 

6.14E+17 

6.09E+1 7 

6.06E+1 7 

6.03E+1  7 

6.03E+1 7 

cc(m/sec)  thermal  speed  of  cathode  gas 

525 

531 

535 

538i 

541 

mdot,cgas(mg/sec) 

5.56E-02 

5.58E-02 

5.60E-02! 

5.59E-02 

5.62E-02 

Vdot  cgas(sccm)  (keep  at  min.=about  2sccm) 

2.00E+00 

2.00E+00 

2.01  E+00 

2.01  E+00 

2.02E+00 

mass  fration  (mdot  cgas/mdotC60) 

1.78E-01 

9.70E-02 

6.32E-02  i 

4.52E-02 

3.46E-02 
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Table  4.2 

13cm  Laboratory  Fullerene  Thruster 


INPUT/CASE  NO 

1 

2 

3 

4 

5 

Dact.(m)  grid  active  dia 

0.13 

0.13 

0.13 

0.13 

0.13 

ds(m)  screen  grid  hole  dia, 

0.002 

0.002 

0.002 

0.002 

0.002 

da(m)  acc.  grid  hole  dia. 

0.002 

0.002 

0.002 

0.002 

0.002 

p(m)  grid  hole  pitch 

0.003 

0.003 

0.003 

0.003 

0.003 

lg(m)  grid  spacing 

0.0020 

0.0020 

0.0020 

0.0020 

0.0020 

Vt(Volt)  total  inter,  grid  voltage 

1000 

1500 

2000 

2500 

3000 

Vn(volts)  net  voltage=mi*c''2/e 

500 

1000 

1500 

2000 

2500 

eatu(-)  propellant  utilization 

0.8 

0.8 

0.8 

0.8 

0.8 

Tn  (k)  guess  at  dis.  cham.  Tvapr-ITERATE 

662 

676 

687 

695 

702 

epsbi(volts)  loss  per  beam  ion 

500 

500 

500 

500 

500 

Vdis. (Volts)  guess  discharge  voltage 

30 

30 

30 

30 

30 

etabase(-),  base  eff. 

0.8 

0.8 

0.8 

0.8 

0.8 

Mw(kg/kgmol)  propellt  mol.  wt. 

720 

720 

720 

720 

720 

dstar  vap(mm)  vaporizer  orifice  dia. 

5.8 

5.8 

5.8 

5.8 

5.8 

Mwcathode(kg/kgmole)  c.  gas  mole  weight 

39.9 

39.9 

39.9 

39.9 

39.9 

fcn=(n  of  Ar)/(n  of  C60)  cathode  flow  fract. 

2.05 

1.13 

0.74 

0.53 

0.41 

(vary  fen  to  get  at  most  3sccm  cathode  flow) 

OUTPUT -THRUSTER 

fo(-)  open  area  fraction 

0.349 

0.349 

0.349 

0.349 

0.349 

Aat(m2)  surface  area  of  grid 

1.327E-02 

1.327E-02 

1.327E-02 

1.327E-02 

1.327E-02 

Aao(m2)  open  area  of  acc.  grid 

4.629E-03 

4.629E-03 

4.629E-03 

4.629E-03 

4.629E-03 

N  (-)numbr  of  holes 

1474 

1474 

1474 

1474 

1474 

mi(kg)  mass  of  prop,  ion 

1.195E-24 

1.195E-24 

1 .1 95E-24 

1 .1 95E-24 

1.195E-24 

ci(m/sec)  ion  exit  speed 

1.158E+04 

1.637E+04 

2.005E+04 

2.31 5E+04 

2.589E+04 

le(m)  effective  grid  spacing 

2.236E-03 

2.236E-03 

2.236E-03 

2.236E-03 

2.236E-03 

J(A/m2)  beam  current  density 

1.289E+01 

2.367E+01 

3.645E+01 

5.093E+01 

6.695E+01 

Ibeam  (A)  beam  current 

0.060 

0.110 

0.169 

0.236 

0.310 

ni(1/m3)  ion  number  density 

6.948E+1 5 

9.025E+1 5 

1.135E+16 

1 -373E+1 6 

1.614E+16 

cn(m/sec)  speed  of  neutrals 

124 

125 

126 

127 

127 

ni/nn(-) 

0.0427 

0.0305 

0.0251 

0.0219 

0.0197 

nn(1/m3)  neutral  number  density 

1.626E+17 

2.957E+1 7 

4.516E+17 

6.275E+1 7 

8.207E+1  7 

ntot(1/m3)  tot  number  density 

1.696E+17 

3.047E+1 7 

4.629E+1 7 

6.412E+17 

8.368E+1 7 

ntot(1/m3)  based  on  Tdis.  chamb.  input 

1.650E+17 

2.940E+1 7 

4.551  E+1 7 

6. 199E+1 7 

8.076E+1  7 

ntot(=R42)/ntot(=R43)lteterate  to  =1 

1.028E+00 

1.036  E+00 

1.017E+00 

1.034E+00 

1.036E+00 

Tdis.  cavity  (degC)  based  ntot 

389 

403 

414 

422 

429 

p(N/m2)  in  dis.  cavity  based  on  ntot 

1.549E-03 

2.843E-03 

4.389E-03 

6.1 50E-03 

8.1 07E-03 

p(Torr)  in  dis.  cavity  based  on  ntot 

1 ,1 54E-05 

2.1 1 8E-05 

3.270E-05 

4.582E-05 

6.040E-05 

F(N)  thrust  from  eq.4 

0.0052 

0.0134 

0.0253 

0.0408 

0.0599 

F(grams))  thrust  from  eq.4 

0.53 

1.37 

2.58 

4.16 

6.11 

mdot(kg/sec)  mass  flow  eq.5 

5.563E-07 

1.022E-06 

1.573E-06 

2.1 99E-06 

2.891  E-06 

mdot(mg/sec)  mass  flow  eq.5 

0.556 

1.022 

1.573 

2.199 

2.891 

Isp(sec)  specific  impls 

947 

1338 

1638 

1891 

2114 

etat(-)  thrust  efficiency 

0.40 

0.53 

0.60 

0.64 

0.67 

Pttl  =  FA2/(2*mdot*eta) 

59.96 

164.98 

338.44 

590.98 

931.98 

Pdis=Pttl-Pbeam=Pttl-lbeam*Vnet 

30.14 

55.41 

85.40 

119.49 

157.24 

ldis.(amps)=Pdis/Vdis 

1.00 

1.85 

2.85 

3.98 

5.24 

OUTPUT -VAPORIZER 

guess  Tv(degC) 

651 

680 

701 

719 

733 

Pc60(mTorr) 

4.07E+01 

7.64E+01 

1.18E+02 

1.68E+02 

2.20E+02 

Pc60(N/m2) 

5.46E+00 

1.03E+01 

1.58E+01 

2.26E+01 

2.95E+01 

mdot  c60  (mg/sec) 

5.57E-01 

1.03E+00 

1.57E+00 

2.22E+00 

2.89E+00 

mdot/mdot(from  row#52) 

(iterate  Tv  to  get  Mdot/mdot=1) 

1.00 

1.01 

1.00 

1.01 

1.00 

OUTPUT  -  cathode  max  flow 

Nc.gas(1/m3)  cathode  gas  number  densty 

3.48E+1 7 

3.44E+1 7 

3.43E+1 7 

3.40E+1 7 

3.43E+1 7 

cc(m/sec)  thermal  speed  of  cathode  gas 

525 

531 

535 

538 

541 

mdot.cgas(mg/sec) 

5.60E-02 

5.60E-02 

5.62E-02 

5.61  E-02 

5.69E-02 

Vdot  cgas(sccm)  (keep  at  min.=about  2sccm) 

2.01  E+00 

2.01  E+00 

2.02E+00 

2.01  E+00 

2.04E+00 

mass  fration  (mdot  cqas/mdotC60) 

1.01  E-01 

5.48E-02 

3.57E-02 

2.55E-02 

1.97E-02 
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The  second  approach,  where  one  can  consider  feeding  pelletized  fuel  using  gun-like 
magazines,  requires  complex  mechanical  system  with  a  hermetic  seal  in  between  the  vaporizer  and 
the  storage  tank  so  that  fullerene  vapor  does  not  flow  back  into  the  storage  tank  where  it  would 
condense.  The  hermetic  seal  would  have  to  be  basically  a  high  temperature  valve  operating  at  up  to 
700°C  so  that  condensation  does  not  occur  within  it,  which  could  freeze  it  up  and  plug  the  conduit. 
The  pros  and  cons  of  these  approaches  are  summarized  in  Table  4.3. 

4.2.1  Laboratory  Vaporizer  and  Storage 

As  indicated  in  Fig.  4.1  for  the  laboratory  version  of  the  FIT-13,  we  selected  Option  #1 
from  Table  4.3  which  is  the  more  likely  approach  for  the  "early"  fullerene  thruster  application  in 
space  and  for  material  processing  on  the  ground. 

Per  Table  4.2  the  design  vaporization  rate  is  5  mg/sec.  The  volume  of  the  vaporizer  was 
selected  to  be  850  cm3  which  will  deliver  5  mg/sec  for  approximately  50  hours  when  initially  60% 
full.  The  maximum  vaporizer  temperature  was  selected  to  be  750°C  because  at  approximately 
800°C  fullerenes  start  to  fragment.  Because  the  fullerene  vapor  pressure  is  exponentially 
dependent  on  temperature,  the  vaporizer  internal  surface  temperature  uniformity  requirement  was 
selected  to  be  ±5°C.  Worst  case  analysis  indicate  that  the  required  heater  power  to  maintain  the 
laboratory  vaporizer  at  the  design  temperature  is  1500  Watts.  This  can  be  mitigated  by  adding 
more  radiation  shields. 

The  fullerene  wetted  inside  surface  of  the  vaporizer  are  made  of  SST  316.  The  temperature 
uniformity  is  achieved  by  explosion  bonding  Cu  layer  into  the  SST  inside  layer.  The  copper  acts 
as  a  heat  spreader.  Heating  was  initially  accomplished  by  tungsten  wire  threaded  through  A  i  2O3 
tubes  and  placed  in  grooves  in  the  Cu  layer.  Later  the  tungsten  wire  and  A^?03  insulators  were 
replaced  by  swaged  tantalum  cable.  The  general  construction  is  shown  in  the  photograph  of  the 
partially  assembled  vaporizer  in  Fig.  4.2. 

Fullerene  vapor  is  transported  through  two  0.375  dia.  SST  304  tubes  into  the  fullerene 
manifold  located  inside  the  discharge  chamber.  The  size  of  the  transport  lines  was  selected  on  the 
basis  of  Fig.  4.3  which  shows  an  orifice  diameter  vs.  fullerene  vapor  temperature.  The  graph  was 
constructed  with  the  aid  of  Eqs.  [3.1]  and  [3.2]. 

4.2.2  Future  Vaporizer  and  Feeding  System  Concept 

The  fundamental  options  of  heating  all  propellant  in  a  single  large  volume  that  serves  as  a 
storage  and  a  vaporizer  or  a  small  fraction  of  the  propellant  fed  incrementally  into  a  small  vaporizer 
were  discussed  at  the  beginning  of  Section  4.2  and  Table  4.3.  For  long  term  missions  it  may  be 
advantageous  to  choose  Option  2  from  Table  4.3  for  the  following  reasons: 

1)  The  C60/C70  thermal  degradation  appears  to  be  a  function  of  not  only  temperature  but 
also  time  at  temperature  and  number  of  reheat  cycles.  Thus,  heating  a  large  volume  of 
propellant  and  keeping  it  hot  for  the  whole  time  an  orbit  may  not  be  advisable  (aside  of 
heat  losses  consideration). 

2)  Heating  a  large  volume  of  propellant  just  for  the  thruster  operating  period  which  may 
be  as  short  as  1/2  hour  per  day  would  be  extremely  energy  intensive  and  may  affect  the 
vapor  pressure  depending  on  the  number  of  reheat  cycles. 

Thus,  feeding  a  small  fraction  of  the  propellant  into  a  small  vaporizer  with  fast  heat  up  and 
low  thermal  losses  may  be  the  best  approach.  Final  determination  will  however,  require  data  from 
long  duration  fullerene  vaporization  experiments.  The  transport  of  the  solid  fullerenes  from  a  large 
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Fig.  4.2  Partially  Assembled  Fullerene  Vaporizer 
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storage  tank  to  a  small  vaporizer  and  the  rejection  of  unvaporized  fullerene  residue  from  the 
vaporizer  which  may  be  necessary,  makes  this  approach  difficult  to  implement.  A  concept  which 
appears  to  deserve  further  exploration  is  shown  in  Fig.  4.4  which  depicts  a  conceptual  sketch  of  an 
incremental  fullerene  storage  and  feeding  system.  The  C(,o  powder  is  stored  pressed  in  a  cartridge 
of  some  desirable  volume,  say  10  to  50  grams  depending  on  the  mission.  An  arbitrary  number  of 
cartridges  are  strung  together  in  a  machine  gun-like  belt.  The  cartridges  could  be  made  of  some 
suitable  light  weight  plastic,  metal  or  even  quartz.  To  insert  the  propellant  into  the  vaporizer,  the 
cartridge  is  positioned  in  front  of  the  breach  and  a  piston  forces  the  compressed  powder  into  a 
small  vaporizer.  The  breach  closes  and  provides  a  gas  tight  seal  at  temperature  T  =  700°C.  If  it  is 
shown  that  unvaporized  residue  is  significant  and  must  be  periodically  removed  from  the 
vaporizer,  then  the  whole  cartridge  (for  example,  made  of  quartz)  must  be  inserted  into  the 
vaporizer  and  heated  and  after  a  cool  down,  the  cartridge  with  the  unvaporized  residue  in  it  is  then 
ejected. 


There  are  many  variations  on  this  basic  concept.  An  appealing  feature  of  this  approach  is 
that  the  unvaporized  residue  (if  any)  can  be  rejected  and  that  the  system  is  scaleable  for  any 
quantity  of  propellant  by  making  the  belt  longer.  The  mechanism  does  not  appear  overly 
complicated  except  for  the  breach  (i.e.,  gas  tight,  high  temperature  valve). 

4.3  Thruster  Design 

The  schematic  of  the  FIT- 13  thruster  is  shown  in  Fig.  4.5.  Unlike  conventional  Xe  fueled 
ion  thruster,  the  FIT  must  have  heated  walls  to  prevent  condensation  of  fullerene  ions.  The 
minimum  wall  temperatures  are  listed  in  Tables  4.1  and  4.2  for  two  different  grid  designs.  Cusp 
magnetic  field  was  initially  provided  by  the  indicated  electromagnets  and  later  replaced  by 
permanent  rod  magnets  (see  Section  5).  Another  consequence  of  high  discharge  chamber  wall 
temperature  is  the  need  to  provide  heat  shields  around  the  whole  vaporize/thruster  flow  train  and  to 
accommodate  differential  thermal  growth  of  the  discharge  chamber  and  the  grids  which  necessitates 
pivoting  dielectric  grid  mounting  system.  These  features  that  set  it  apart  from  the  conventional  ion 
thruster  and  the  salient  components  of  the  FIT-13  are  described  in  the  subsequent  section. 

4.3.1  Cathode 

The  hollow  cathode  schematic  is  shown  in  Fig.  4.6.  The  fundamental  approach  to  the 
design  follows  NASA  JPL  and  NASA  Lewis  designs. (21>22)  The  cartridge  emitter  is  tungsten 
impregnated  with  Ba  and  Ca  oxides.  The  orifice  size  is  approximately  0.033".  The  orifice  plate  is 
made  of  W  and  replaceable  to  allow  easy  change  of  orifice  size.  The  cathode  shroud  (enclosure)  is 
isolated  and  serves  as  a  keeper  electrode.  The  heater  was  initially  made  of  BN  with  VM  type 
tungsten  wire  as  the  heating  element  and  later  replaced  by  swaged  tantalum  cable.  This  approach 
was  taken  to  permit  the  replaceability  of  the  orifice  plate  size.  Orifice  size  in  combination  with 
internal  Ar  pressure  dictates  the  Ar  mass  flow  which  must  be  sufficient  to  preclude  diffusion  of 
C6o  ions  into  the  hollow  cathode  interior  thereby  lowering  its  emissivity. 

Several  tests  were  conducted  with  the  hollow  cathode  replaced  by  2%  Th  tungsten  wire 
cathode  that  was  heated  by  separate  power  supply.  The  design  of  this  cathode  is  shown  in 
Fig.  4.7. 

4.3.2  Discharge  Chamber 

The  unconventional  features  of  the  FIT- 13  are  the  discharge  chamber  heaters  and  the 
magnetics  to  create  the  desired  cusp  field. 
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The  heaters  initially  consisted  of  a  double  wound  coil  made  of  VM  tungsten  wire  with 
radiation  heat  shields.  There  were  two  identical  heaters  that  fit  over  the  outside  diameter  of  the 
discharge  chamber  as  indicated  in  Fig.  4.5.  The  double  winding  was  employed  to  ensure  that 
magnetic  fields  generated  by  the  heaters  are  near  zero  as  one  winding  cancels  the  magnetic  field  of 
the  other.  The  estimated  total  heat  power  to  maintain  the  discharge  chamber  walls  at  about  500°C  is 
350  Watts.  This  assumes  no  discharge  on  the  inside  and  which  reduces  the  heater  input  power. 
Further  reduction  can  be  obtained  by  additional  heat  shields.  The  VM  tungsten  heaters  were  later 
replaced  by  a  heated  insert  into  the  discharge  chamber  which  is  further  discussed  in  Section  5. 

The  magnetic  structure  is  shown  in  Fig.  4.8.  The  goal  of  the  design  was  to  create  magnetic 
field  contours  similar  to  the  high  performance  thrusters  with  cusp  field  reported  by  Sovey(20)  and 
Matossian,  et  al.C21)  The  field  was  initially  created  by  a  set  of  12  coils  located  on  12  soft  iron  rods. 
The  rods  mate  up  against  a  soft  iron  mid  flange  conducting  the  opposite  direction  flux  into  this 
flange.  The  flux  lines  split  into  B1  and  B2  components  on  the  inside  of  the  discharge  chamber  as 
indicated.  The  B1  flux  lines  are  conducted  through  the  soft  iron  flange  that  surrounds  the  hollow 
cathode  in  the  rear  of  the  thruster.  The  B2  flux  lines  are  conducted  through  the  front  soft  iron 
flange.  Other  magnetic  materials  may  be  used  in  the  flight  version  of  the  thruster  that  have  higher 
saturation  field  and  therefore  require  less  material  to  conduct  the  flux.  A  list  of  materials  to  be 
considered  is  shown  in  Table  4.4.  The  electromagnetic  coils  were  later  replaced  by  permanent 
magnets  which  are  further  discussed  in  Section  5. 

The  design  calculations  were  carried  out  using  the  2D  version  of  a  commercial  finite 
element  numerical  code  known  as  Maxwell.  The  geometry  of  the  physical  situation  which  as 
implemented  is  three-dimensional,  had  to  be  modified  for  the  2D  software  package.  The  modified 
geometry  which  approximates  the  actual  design  in  the  thruster  interior  flux  and  B  field  magnitude  is 
shown  in  Fig.  4.9.  The  12  magnets  employed  in  practice  are  replaced  by  2  coils  as  indicated  and 
the  magnetic  rods  are  replaced  by  a  continuous  soft  iron  shell.  The  cross-sectional  area  of  each 
physical  rod  is  equal  to  1/6  of  the  cross-sectional  area  of  the  shell.  Matching  the  flux  in  the  shell  to 
the  iron  rods  will  yield  the  same  B  field  distribution  in  the  discharge  chamber.  The  calculated  flux 
lines  are  shown  in  Fig.  4.10  and  the  maximum  field  contours  are  shown  in  Fig.  4.11. 

4.3.3  Grid  Design 

Two  grid  sets  have  been  designed  following  the  principals  set  out  by  Rovang  and 
Wilbur.(22)  Each  set  is  specified  in  Table  4.5.  The  major  difference  between  the  two  sets  is  the 
hole  pitch  (center-to-center  distance)  which  is  4  mm  for  set  1  and  3  mm  for  set  2.  As  a  result  set  1 
has  829  holes  in  each  grid  and  set  2  has  1474  holes  in  each  grid.  The  material  of  construction  is 
Molybdenum.  To  save  funds  the  grids  were  made  flat  and  only  one  set  has  been  made. 
Performance  of  each  set  is  estimated  in  Tables  4.1  and  4.2  and  also  in  Appendix  A. 

Each  grid  is  mounted  on  6  pivoting  posts  made  of  hot  pressed  SiN.  Each  post  pivots  on  a 
pin  with  its  retaining  bracket  mounted  on  the  discharge  chamber.  This  allows  unrestrained 
independent  differential  thermal  expansion  in  the  radial  direction  of  the  discharge  chamber  and  the 
grids. 


4.3.4  Electrical  System 

The  FIT-13  electrical  system  schematic  is  shown  in  Fig.  4.12.  It  appears  unusually 
complex  because  it  contains  several  elements  not  associated  with  a  conventional  ion  thruster.  This 
includes  the  fullerene  vaporization  chamber  with  two  heaters  each  with  an  independent  power 
supply,  discharge  cavity  heater  and  a  power  supply  plus  discharge  chamber  electromagnets  with  its 
power  supply.  Additionally,  the  schematic  contains  separately  heated  beam  target/collector. 
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Fig.  4.8  Magnetic  Structure  for  the  FIT- 13  Laboratory  Thruster 
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Table  4 -5  FIT-13  Grid  Specifications 


Parameter 

Set  1 

Set  2 

da  (mm) 

2 

2 

ds  (mm) 

2 

2 

p  (mm) 

4 

3 

lg  (mm) 

1 

1 

ts  (mm) 

0.5 

0.5 

ta  (mm) 

0.5 

0.5 

N  (number  of  holes) 

829 

1474 

Pattern 

Rectangular 

Rectangular 

Active  Diameter  (mm) 

130 

130 
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All  power  supplies  are  floating  fed  with  separate  isolation  transformer.  The  entire  system 
is  remotely  controlled  from  a  PC  style  computer  which  also  functions  as  the  data  acquisition 
system.  Further  description  of  this  system  is  relegated  to  Section  6. 
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5.0  THRUSTER  CONSTRUCTION  AND  DESIGN  EVOLUTION 


All  machined  components  of  the  thruster  system  were  fabricated  by  outside  vendors  per 
Busek  drawings  and  assembled  in-house.  Photographs  of  the  major  components  are  shown  in 
Figs.  5.1  through  5.3  showing  the  vaporization  chamber,  the  back  flange  of  the  thruster  and  the 
discharge  chamber,  respectively.  The  vaporizer  design  and  construction  features  were  discussed  in 
Section  4.  A  key  aspect  is  internal  temperature  uniformity  requirement  achieved  through  the  use  of 
explosion  bonded  Cu  and  SST  vaporizer  walls.  The  back  flange  of  the  thruster  (Fig.  5.2)  is  made 
of  soft  iron  acting  as  a  guide  of  the  magnetic  field  flux.  The  flange  also  contains  the  fullerene 
distribution  manifold  made  of  SST.  The  photograph  of  the  discharge  chamber  (Fig.  5.3)  shows 
the  mid  and  front  magnetic  flanges  to  form  the  internal  cusp  magnetic  field  structure.  Also  shown 
is  a  portion  of  the  chamber  wall  heaters  to  prevent  internal  fullerene  condensation  and  pivoting  grid 
mounting  system.  The  pivoting  insulators  made  of  SiN  which  allows  free  radial  and  azimuthal 
differential  thermal  expansion  of  both  the  chamber  and  the  grids. 

Photographs  of  the  assembled  thruster  system  are  shown  in  Figs.  5.4  and  5.5.  The 
connection  between  the  vaporizer  and  the  discharge  chamber  is  revealed  in  Fig.  5.4  along  with  a 
portion  of  the  hollow  cathode  mounted  on  a  BN  flange.  The  electromagnet  coils  creating  the 
internal  cusp  field  are  visible  in  both  Figs.  5.4  and  5.5,  the  later  also  shows  the  grid  and  the  grid 
mounting  system.  The  entire  assembly  is  mounted  on  ceramic  standoffs  which  in  turn  are  mounted 
on  a  grounded  skid  with  rails.  The  rails  support  grounded  shroud  (not  shown)  that  encloses  the 
entire  assembly  and  acts  as  the  outer  most  heat  shield. 

Following  assembly,  the  magnetic  field  prediction  (see  Section  4)  was  verified  using  a 
Gauss  meter  and  the  traditional  technique  of  flux  mapping  using  iron  filings.  Figure  5.6  shows  a 
photograph  of  the  flux  lines  formed  by  iron  filings  on  cardboard  paper  inserted  into  the  discharge 
chamber.  This  photograph  and  measurements  confirm  predictions  shown  in  Figs.  4.10  and  4.11. 

Initial  experiments  revealed  several  difficulties  with  operating  the  thruster  and  led  to  the 
following  design  changes. 

1)  The  external  discharge  chamber  heater  was  replaced  by  heated  sleeve  that  is  inserted 
into  the  discharge  chamber.  The  design  of  this  sleeve  is  shown  in  Fig.  5.7.  Heating  is 
accomplished  by  coaxial  swaged  wire  whose  outer  sheet  is  periodically  tack  welded  to 
the  sleeve. 

2)  The  change  in  item  1  facilitated  improved  discharge  chamber  heating  and  permitted 
replacement  of  electromagnets  Alnico  5,  1/2"  dia.  cylindrical  permanent  magnets  that 
yielded  approximately  the  same  magnetic  field. 

3)  The  stainless  steel  conduits  for  fullerene  vapor  that  connect  the  vaporizer  and  the 
discharge  chamber  were  replaced  with  copper  tubes  of  the  same  size  to  affect  better 
temperature  uniformity  and  prevent  internal  C<30  vapor  condensation. 

4)  Initially  tested  hollow  cathode  heater  was  later  replaced  by  swaged  tantalum  cable. 

5)  Initially  tested  tungsten  wire  vaporizer  heater  was  replaced  by  swaged  tantalum  cable. 

The  results  reported  in  subsequent  sections  were  obtained  after  implementation  of  the  above 
discussed  changes. 
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Fig .  5.4  Side  view  of  the  FIT-13  thruster  system  showing  C^o  vaporization  (left) 
and  the  discharge  chamber  (right) 
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Fig.  5.6  Iron  Fillings  Indicate  Distribution  of  the  Field  Lines  Inside  the  Discharge  Chamber 


0266 


45 


BUSEK 


BUSEK  Report  No.  036-3 


6.0  FACILITY 


The  vacuum  tank  is  a  cylindrical  stainless  steel  vessel  approximately  4  ft.  in  diameter  and  5 
ft.  long  with  hinged  doors  on  both  ends  of  the  cylinder  and  3  flanged  access  ports  each  24"  dia. 
The  tank  is  pumped  with  a  24"  oil  diffusion  pump  backed  by  355  CFM  two-stage  system  that 
consists  of  a  blower  and  mechanical  pump. 

The  thruster  system  is  installed  within  a  20"  dia.  cylindrical  side  enclosure  as  shown  in 
Fig.  6.1.  The  thruster  system  is  mounted  on  a  grounded  bracket  within  this  enclosure  and 
electrically  floats  at  the  beam  power  supply  potential.  In  order  to  prevent  streaming  of  neutralizer 
electrons  to  the  thruster  system  a  grounded  baffle  is  installed  near  the  acceleration  grid  as  indicated 
in  Fig.  6.1.  A  commercial  neutralizer  cathode  is  operated  on  Ar  as  is  the  internal/discharge  cathode 
when  in  use. 

The  electrical  system  consists  of  9  dc  power  supplies,  installed  in  3  adjacent  24" 
instrumentation  cabinets.  Each  power  supply  can  float  with  respect  to  ground  and  can  be 
connected  to  arbitrary  potential  in  the  system  (see  Fig.  4.11).  Isolation  rated  to  3  kV  is 
accomplished  by  feeding  each  power  supply  with  an  isolation  transformer  and  mounting  the  power 
supply  on  fiberglass  rails  inside  the  grounded  metallic  cabinets.  An  overall  photograph  of  the 
system  is  shown  in  Fig.  6.2. 

Control  of  the  power  supplies  and  data  acquisition  is  performed  by  a  PC  style  computer. 
The  interface  between  the  floating  high  voltage  system  and  the  grounded  computer  consists  of  high 
voltage  (3  and  4.5  kV)  isolation  amplifiers  mounted  on  two  different  types  of  printed  circuit 
boards.  One  type  is  used  for  straight  isolation  without  amplification  (Datatrak  boards,  4.5  kV 
isolation)  and  the  other  type  is  used  for  thermocouple  signals  and  has  adjustable  gain  (KD  Physics, 
3  kV  isolation).  The  system  is  controlled  and  data  are  stored  using  software  called  LabView  made 
by  National  Instruments.  Operator  inputs  his  commands  either  via  the  PC  keyboard  or  a  mouse. 
Each  power  supply  can  also  be  controlled  manually  by  using  isolated  control  knob  extension  shafts 
to  ensure  operator  safety.  The  manual  control  and  manual  data  taking  was  often  used  instead  of  the 
computerized  system  because  of  frequent  failures  of  the  isolation  amplifiers  that  shorted  when  the 
grids  sparked. 
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Fig.  6.2  Large  Throughput  Vacuum  Facility 
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7.0  EXPERIMETNAL  RESULTS  &  DISCUSSION 


7.1  Procedures 

Fullerene  propellant  was  prepared  two  different  ways.  The  simpler  preparation  method 
involved  simple  vacuum  drying  of  the  as  purchased  sample  while  the  more  complex  process 
involved  sublimation  of  the  as  purchased  sample,  capture  and  condensation  of  the  vapor.  Only  the 
condensate  was  then  used  as  the  propellant.  The  propellant  typically  consisted  of  C60/C70  mix 
with  80/20  split.  The  suppliers  and  preparation  details  were  discussed  in  detail  in  Section  3.0. 

Regardless  of  the  preparation  procedure  the  propellant  was  always  weighed  and  placed  into 
the  thruster  system  vaporizer  prior  to  the  installation  of  the  system  in  the  vacuum  chamber  (see 
Section  6,  Fig.  6.1).  Typical  propellant  batch  had  10  grams.  Upon  test  completion  the  mass  of 
the  propellant  remaining  in  the  vaporizer  was  weighed  were  any  remnants  in  the  discharge 
chamber.  Following  thruster  system  installation  into  the  test  chamber,  the  chamber  was  pumped 
down  to  low  10'5/10'6  torr  level.  Argon  was  then  introduced  into  the  discharge  chamber  via  the 
hollow  cathode.  Following  cathode  insert  heating  the  cathode  was  started  using  the  keeper  circuit 
followed  by  discharge  initiation  on  the  cathode  Ar  flow  only.  The  Ar  flow  rate  was  generally  kept 
to  a  minimum  as  dictated  by  the  cathode  operation.  Ar  beam  was  established  by  applying  voltages 
to  the  screen  and  the  acceleration  grid.  Typical  screen  voltage  also  referred  to  as  the  beam  voltage 
was  1  kV  (relative  to  ground)  while  the  accelerator  voltage  was  typically  at  -500  V.  Beam 
neutralization  was  achieved  using  a  commercially  available  cathode  operated  on  argon. 

Once  the  thruster  system  was  operating  on  argon  and  the  baseline  Ar  beam  was  established, 
the  discharge  chamber  heating  is  applied  followed  by  heating  of  the  vaporizer.  As  the  fullerene 
mass  flow  increased  the  beam  gradually  changed  appearance  from  blue/purple  to  a  light  yellow 
which  is  a  visual  indication  that  the  discharge  operates  on  fullerenes.  This  is  shown  in  Fig.  7.1 
where  the  top  photo  shows  the  discharge  operating  on  Ar  (as  viewed  through  the  grids)  and  the 
bottom  photo  shows  the  discharge  color  when  operating  on  fullerenes.  The  bright  region  in  the 
center  is  the  dense  Ar  plasma  emanating  from  the  hollow  cathode. 

Several  tests  were  conducted  with  the  Ar  hollow  cathode  replaced  by  2%  Th  heated 
tungsten  filament.  During  these  tests  the  procedure  described  above  was  slightly  modified. 

Several  tests  were  also  conducted  with  stainless  steel  sheet  metal  target  exposed  to  the  beam 
to  explore  carbon  film  deposition.  The  target  was  approximately  30  cm  downstream  from  the  grids 
and  was  always  kept  at  zero  potential  (ground  level).  Results  of  such  test  are  discussed  in  Section 
8.4. 

7.2  Representative  Experimental  Data 

Typical  data  sheets  that  contain  all  measured  quantities  recorded  as  a  function  of  time  are 
shown  in  Tables  7.1  and  7.2.  The  tank  pressure,  cathode  flow,  keeper  current  and  voltage,  the 
neutralizer  flow,  keeper  current  and  voltage,  the  discharge  current  and  voltage,  the  potential  of  the 
grids,  the  accelerator  current  and  the  beam  current  were  measured.  Additionally  the  outputs  of 
three  thermocouples  were  recorded  as  were  the  voltages  and  currents  of  the  vaporizer  and 
discharge  chamber  heaters.  The  thermocouple  labeled  Tel  is  located  on  the  front  flange  of  the 
thruster  just  upstream  of  the  grids  (external  to  the  discharge  chamber)  and  is  the  coldest  portion  of 
the  system.  Thermocouple  Tc4  is  located  on  the  vaporizer  manifold  between  the  vaporizer  and  the 
thruster  and  approximately  corresponding  to  the  temperature  of  the  vapor  conduits  between  the 
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Fig.  7.1  Photo  of  the  Discharge  Through  the  Grids  Showing  (a)  Argon  Plasma  Only 
(b)  Fullerene/ Argon  Plasma.  The  former  is  blue/purple,  the  latter  is  yellow. 
The  bright  spot  in  the  center  is  the  dense  Ar  plasma  emanating  from  the  hollow 
cathode 
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(I  e  at  following  Beam  Voltages  1.000V,  800V.  650V.  500V) 
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Few  unsuccesful!  runs  attempted  betore  this  one 

Test  Conditions.  1)  Intergrid  distance  increased  to  0  133"  2)  Disch  Cavity  Heater  Polarity  Changed.  3)  tOg  of  new  Hoechst  batch  C60  loaded  on  6/6/96 
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vaporizer  and  the  thruster.  The  thermocouple  Tc6  is  on  the  inside  of  the  vaporizer.  Tc6  was  used 
to  calculate  the  fullerene  vapor  pressure  (not  measured  directly)  and  consequently  the  mass  flow 
through  the  propellant  orifice  (see  Fig.  4.1  for  orifice  location).  Flow  through  the  orifice  lowered 
the  vapor  pressure  and  Tc4  was  maintained  at  or  above  condensation  temperature  at  that  pressure. 
Tel  was  maintained  at  or  above  the  fullerene  condensation  temperature  that  corresponds  to  the 
anticipated  fullerene  vapor  pressure  in  the  discharge  chamber. 

7.3  Typical  Post  Test  Observations 

The  walls  of  the  discharge  chamber  following  a  typical  test  were  covered  with  thin,  hard 
carbon  film.  This  film  adhered  well  to  the  walls  until  after  several  tests  it  accumulated  to 
measurable  thickness  (~0.1  mm  or  less)  at  which  point  it  formed  poorly  adhering  flakes  that  could 
be  scraped  off  with  a  hard  object.  The  flakes  were  insoluble  in  benzene  and  appeared  to  be 
consistent  with  diamond-like  carbon  (DLC).  There  was  no  film  in  the  locations  of  highest 
magnetic  flux  where  the  discharge  current  electrons  are  likely  to  strike  the  chamber  walls  (anode). 
The  film  free  areas  were  bear  metal  bands  around  the  circumference  of  the  discharge  chamber  about 
6  mm  wide  in  the  locations  of  the  external  magnetic  flanges. 

The  upstream  side  of  the  screen  grid  was  coated  with  extremely  hard  carbon  film  which 
could  be  removed  with  difficulty  only  by  vigorous  sanding.  The  film  appeared  to  be  hardest  and 
best  adhering  in  the  center  of  the  screen  grid  covering  an  area  about  30  to  40  mm  in  diameter. 
Wetting  the  film  with  benzene  did  not  dissolve  it.  It  also  appeared  as  a  DLC  film  although  no 
quantitative  analysis  were  performed. 

This  observation  is  consistent  with  DLC  film  deposition  by  Gruen(-3)  who  used 
argon/fullerene  plasma  created  by  microwave  discharge. 

The  upstream  side  of  the  acceleration  grid  was  coated  with  thin  brownish  film  primarily  in 
the  areas  where  the  screen  and  accelerator  grid  misaligned  due  to  their  differential  expansion  (see 
Fig.  7.1).  This  film  was  soft  and  could  be  wiped  off. 

For  tests  utilizing  the  hollow  cathode,  there  was  no  significant  carbon  debris  in  the  chamber 
however,  after  test  with  the  filament  cathode  (both  with  and  without  Ar)  there  was  copious  amount 
of  carbon  in  the  chamber  approaching  the  total  vaporized  fullerene  mass.  This  is  consistent  with 
the  Anderson's,  et  al(7)  results  who  reported  up  to  2/3  of  the  vaporized  fullerenes  remaining  in  the 
discharge  chamber  as  fullerene  fragments  when  using  filament  cathode. 
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8.0  DATA  ANALYSIS  &  DISCUSSION 


The  equations  used  to  reduce  the  experimental  data  are  summarized  in  Table  8.1.  All  of 
these  equations  are  straightforward,  conventional,  relationships  that  require  no  explanations  except 
for  two  aspects  related  to  the  estimation  of  the  average  intergrid  gap  and  the  appropriate  ion  mass 
that  is  required  for  evaluation  of  the  space  charge  limited  beam  current  density  and  thrust.  These 
two  issues  are  discussed  next  followed  by  the  presentation  of  the  experimental  results. 

8.1  Average  Intergrid  Gap 

The  intergrid  separation,  is  not  precisely  known  due  to  thermal  distortion  of  the  grids.  For 
the  test  in  Table  7.1  the  intergrid  distance  was  estimated  to  be  1.91  mm  (see  bottom  right  corner  of 
the  table).  As  indicated  in  Section  4,  the  grids,  due  to  high  cost  of  dished  grids,  were  flat  0.5  mm 
thick  Molybdenum  and  allowed  to  grow  radially  on  pivoted  insulators.  Because  the  screen  grid  is 
at  higher  temperature  than  the  accelerator  (outer)  grid  it  bows  outward  more  than  the  accelerator 
grid  as  schematically  indicated  in  Fig.  8.1. 

During  our  initial  tests,  the  cold  state  intergrid  distance  was  set  to  2  mm  which  always 
resulted  in  grid  electrical  shorting  while  approaching  the  operating  temperature.  This  led  us  to 
increase  the  cold  separation  to  between  2.70  to  3.38  mm.  At  2.7  mm  the  grids  still  occasionally 
shorted  so  that  the  minimum  hot  condition  intergrid  gap  at  the  center  was  estimated  to  be  3.38  - 
2.70  =  0.62  mm.  The  average  gap  shown  in  Tables  7.1  and  7.2  was  then  estimated  using  the 
following  weighted  average  gap  relationship 
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8.2  Beam  Ion  Mass 

The  space  charge  limited  beam  current  density  (JBCL>  Eq.  18  in  Table  8.1)  and  the  thrust 
(T,  Eq.  20  in  Table  8.1)  depend  on  ion  mass  (mi).  When  the  discharge  operates  on  a  mixture  of 
Ar,  C50  and  C70  it  is  not  immediately  evident  which  ion  species  will  dominate  the  beam.  To 
simplify  the  data  reduction  we  discuss  the  limiting  cases,  i.e.,  when  the  beam  consists  of  argon 
ions  and  when  it  consists  of  predominantly  fullerene  ions. 

Before  the  vaporizer  is  heated  to  above  (~400°C)  the  C^q  +  C70  vapor  flow  is  negligible 
and  the  discharge  operates  on  Ar  (supplied  via  the  hollow  cathode  or  in  the  case  of  filament  cathode 
via  an  adjacent  port).  Thus  in  this  case  the  appropriate  mj  is  that  corresponding  to  Ar  (mjAr  =  6.63 

x  10  26  kg).  When  the  vaporizer  is  delivering  mass  flow  of  fullerenes  (mf)  that  is  comparable  or 

greater  than  the  flow  of  argon  mass  flow  (mAr),  the  ion  number  density  of  each  species  in  the 
discharge  chamber  and  the  beam  composition  is  dominated  by  fullerene  ions.  This  is  supported  by 
the  following  analysis. 

Consider  the  classical  experiments  carried  by  Rapp,  et  al.,(24)  where  they  measure  the  total 
ionization  cross-sections  (Q)  by  measuring  electron  beam  current  (Ie)  that  intersects  a  molecular  (or 
atomic)  beam  and  produces  an  ion  beam  (1)).  The  cross-section  is  then  computed  as 
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Fig.  8.1  Schematic  of  the  Thermal  Distortion  of  the  Grids 
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TABLE  8.1  (CON'T) 
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e  =  — — 

n  L 


where  n  is  the  neutral  gas  number  density  and  L  is  the  effective  path  length  of  the  electrons  through 
the  neutral  beam  which  impact  ionize  the  neutrals.  Conceptually  Rapp's(24)  experiment  is  the  same 
as  our  discharge  chamber  that  contains  both  Ar  and  Cfto-  Therefore,  the  respective  ion  currents  can 
be  estimated  using  Eq.  [8.1]  as 

4  60  ^e(Qn^)c60  Qcbi)\(  rooi 


4 ,  le.JUJ  L 

Within  the  discharge  chamber  but  away  from  the  electrodes  and  away  from  sheath  region  of 
the  grids  the  electrical  fields  are  low  and  the  ions  will  drift  at  essentially  thermal  velocity. 
Therefore  the  ion  current  ratio  from  Eq.  [8.2]  can  be  also  written  as 


(«m4 
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(TT4 


where  uj  is  the  thermal  velocity  of  the  respective  ion  ^ T /  mi ) .  Within  the  sheath  the  ions  will 

acquire  Bohm  velocity  which  however,  has  the  same  dependence  on  mj.  Assuming  the  same 
temperature  of  both  types  of  species  and  equating  Eqs.  [8.2]  and  [8.3]  yields  the  ratio  of  the  ion 


number  densities 
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where  mc60  and  mAr  are  weights  of  the  C60  molecule  and  argon  atom,  respectively. 

To  estimate  the  neutral  number  densities  ncoo/nAr  we  use  perfect  gas  law  p  =  nkT  and 
known  mass  flows  ( m ) 

n,JL,±*fML  [8.5] 

"  kT  kT  A 'll  M„  11 

where  p  is  the  partial  pressure  of  Ar  or  C60,  A  is  some  flow  area,  Mw  is  molecular  weight  and  k 
and  R  are  Boltzman  and  universal  gas  constants,  respectively.  Again,  assuming  that  the  Ar  and 
C60  neutrals  are  at  the  same  temperature,  the  ratio  of  neutral  number  densities  can  be  estimated 
from  Eq.  [8.3]  as 


nAr  ,  mAr  {  M»C60 


The  ion  currents  in  Eq.  [8.3]  can  be  viewed  as  the  flux  of  charges  toward  the  screen  grid 
where  they  "feel"  the  applied  electric  field  and  get  extracted.  This  flux  can  be  estimated  in  terms  of 
known  quantities  by  combining  Eqs.  [8.2]  and  [8.6]  as 
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By  combining  Eqs.  [8.4]  and  [8.6]  one  obtains  the  ion  densities  ratio 


[8.7] 
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Therefore,  the  ratio  of  the  ion  number  densities  is  only  dependent  on  the  mass  flow  of  the 
respective  gases  entering  the  discharge  chamber  and  the  respective  ionization  cross-sections.  The 
flux  of  the  charges  toward  the  sheath  grid  which  approximates  the  extractable  beam  currents  are 
obtained  by  combining  Eqs.  [8.7]  and  [8.8] 
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which  depends  only  the  ion  densities  and  the  respective  weights. 

To  assess  the  value  of  ion  density  ratio  and  extractable  beam  currents,  the  ionization  cross- 
sections  of  Ar  and  C60  were  taken  from  Refs.  24  and  25  at  electron  energy  of  38  eV 
(approximating  discharge  voltage)  yielding 

Q»  -  233xl0"6cm2  _om3 

Qc60  53. 5 xlO"16  cm2 


Thus,  according  to  Eqs.  [8.8]  and  [8.9]  with  typical  mass  flow  ratio  mAr  /  mC60  <  1,  the  argon  ion 
density  is  negligible  and  the  extractable  argon  beam  is  also  small  relative  to  the  C<$o  beam. 

This  argument  is  further  reinforced  by  considering  the  Penning  ionization  reaction  and 
charge  transfer  discussed  in  Section  2.5  which  indicates,  that  argon  ion  (Ar+)  is  short  lived  within 
the  C60  and  C70  environment  provided  that  Ar+  mean  free  path  for  Ar+/C6o  collision  is  smaller 
than  the  characteristic  dimension  of  the  discharge  chamber. 

Based  on  the  above  discussion  the  criteria  somewhat  arbitrarily  selected  for  the  ion  mass 
(mi  or  <mi>)  value  to  be  substituted  in  the  equation  of  space  charge  current  density  limit,  thrust 
and  mass  flow  of  ions  (Eqs.  [18],  [20]  and  [22]  in  Table  8.1)  is 

>  2.5  then  mi  =  (/«■}  of  fullerenes  given  by  Eq.  9  in  Table  8.1 

[8.10] 

<  2.5  then  mi  =  miKr 
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8.3  Data  Analysis 

Data  from  3  representative  tests  were  analyzed.  The  first  two  tests  we  will  discuss  used 
argon  fed  hollow  cathode  while  the  last  used  heated  tungsten  filament  cathode.  Before  presenting 
the  analysis  we  describe  using  the  Test  #061796  results  as  an  example  how  the  raw  data  were 
processed. 

The  direct  measurements  given  in  Table  7.1  are  shown  graphically  in  Fig.  8.2.  Parameters 
of  interest  obtained  from  the  direct  measurements  are  shown  in  Fig.  8.3.  All  variables  in  Figs.  8.2 
and  8.3  are  plotted  against  time  (in  minutes)  at  intervals  when  data  were  taken.  This  facilitates  easy 
observation  of  the  impact  of  the  test  engineers  actions  such  as  adjustments  in  discharge  current  or  a 
grid  voltage  on  performance  parameters.  All  variables  are  defined  in  Table  8.1.  Below  we 
describe  how  the  graphs  in  Fig.  8.3  were  prepared  which  applies  to  all  other  tests  to  be  presented 
later. 


Constant  Discharge  Current  Test  (#061796) 

The  fullerene  mass  flow  ( m,  =  mdotf)  using  Eqs.  1  through  8  in  Table  8.1  is  shown  as 
Fig.  8.3a  along  with  the  cumulative  fullerene  mass  obtained  by  integration  of  mf.  The  ratios  of 
mf  / mAr  and  mi/mf  is  plotted  in  Fig.  8.3b.  The  first  ratio  is  used  as  a  criteria  for  calculating  the 

appropriate  value  for  the  beam  ion  mass  (see  Eq.  [8.10])  and  the  second  ratio  is  equivalent  to 
propellant  utilization  using  the  fullerene  flow  rates  only.  Figure  8.3c  compares  the  measured  beam 
current  (Ib),  which  is  the  same  as  that  shown  in  Fig.  8.3e,  to  the  theoretical  beam  current  (IBCL 
calculated  from  the  space  charge  current  limit  (Eqs.  18  and  19  in  Table  8.1).  Figure  8.3d  shows 
the  beam  ion  energy  cost  (eg)  calculated  from  Eq.  14  in  Table  8.1.  Figure  8.3e  shows  the 
propellant  utilization  efficiency  (r]u)  and  the  electrical  efficiency  Cqe)  calculated  using  Eqs.  14  and 
16  in  Table  8.1.  Both  equations  used  the  measured  beam  current,  discharge  power  and  beam 
power  and  use  calculated  ion  mass  and  total  mass  flow  (Eqs.  9  and  11  in  Table  8.1).  Figure  8.3f 
shows  thrust  (T)  and  ISp  calculated  using  Eqs.  20  and  21  in  Table  8.1  Note  that  the  thrust  is 
calculated  using  measured  quantities  and  excludes  the  thrust  contribution  from  flow  of  neutral 
species. 

Having  described  how  each  graph  was  prepared  let  us  now  discuss  the  results.  As  seen  in 
Fig.  8.2a  the  discharge  current  was  held  constant  for  the  duration  of  the  test  as  were  the  beam 
voltage  (Vbeam  =  voltage  on  the  screen  grid  with  respect  to  ground)  and  the  accelerator  voltage 
(Vaccel)  in  Fig.  8.2c.  The  operator  controlled  parameters  that  were  changed  are  the  cathode  Ar 
flow  (Fig.  8.2b)  and  the  vaporizer  temperature  (TVap)  in  Fig.  8.2f.  As  seen  from  Fig.  8.2a  the 
changes  in  the  cathode  flow  (the  only  source  of  Ar)  resulted  in  slight  variation  of  the  discharge 
voltage  and  large  variations  in  the  accelerator  impingement  current  (Iaccel)  shown  in  Fig.  8.2e. 
The  beam  current  slightly  increases  when  the  Ar  cathode  flow  decreases  at  67  minutes  into  the  test 
and  starts  decreasing  when  significant  amount  of  fullerenes  flow  into  the  discharge  (time  >110 
min).  A  decrease  in  Ib  is  expected  reflecting  the  change  over  from  Ar  beam  to  the  heavier  mass 
fullerene  beam.  Figure  8.3a  shows  the  fullerene  flow  rate  (mdotf)  and  its  integral  giving  the 
cumulative  fullerene  flow  (mcummf)  which  indicates  that  10  grams  of  fullerenes  were  vaporized. 
Since  10  gram  wire  initially  loaded  into  the  vaporizer  the  calculations  agree  well  with  the  post  test 
measurement  which  resulted  in  9.86  grams  evaporated  with  0.14  grams  unvaporized  residue. 
Figure  8.3b  shows  the  ratio  of  the  fullerene  to  argon  mass  flow  and  the  ratio  of  the  mass  flow  of 
ions  to  masss  flow  of  fullerenes  The  beam  should  already  be  composed  mostly  out  of  fullerenes 
when  mf  (in  Fig.  8.3a)  reaches  0.4  mg/sec  about  120  minutes  into  the  test.  At  that  point  is 
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Fig.  8.2  Discharge  Current,  Voltage,  Argon  Mass  Flow,  Grid  Potentials,  Beam  Current, 
Accelerator  Grid  Current  and  the  Vaporizer  Temperature  were  Measured  and 
Plotted  as  a  Function  of  Time  into  Test  (in  minutes)  for  Test  ID  #061796 
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Fig.  8.3  Fullerene  Mass  Flow  (mdotf),  Cumulative  Fullerene  Mass  Vaporized  (mcummf), 
Ratio  of  Fullerene  Mass  Flow  to  Argon  Mass  Flow  (mdotf/mdotAr),  Ratio  of 
Fullerene  Ion  Mass  Flow  to  Fullerene  Mass  Flow  (mdoti/mdotf),  the  as 
Measured  Ion  Beam  Current  (Ib),  the  Theoretical  Space  Charge  Limited  Beam 
Current  (Ibcl)>  the  Beam  Ion  Cost  (sb),  the  Electrical  (etae)  and  Propellant 
Utilization  Efficiency  (etau),  and  Finally  the  ISp  and  Thrust  are  Plotted  as  a 
Function  of  Test  Time  in  Minutes.  The  above  parameters  were  calculated  from 
the  measured  data  in  Fig.  8.2  using  Equations  in  Table  8.1. 
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m.J  mf  about  0.5  indicating  that  approximately  half  of  the  fullerene  molecules  are  extracted  as  ion 
beam  current. 

Figure  8.3c  compares  the  theoretical  space  charge  limited  beam  current  (IBCL)  to  the 
measured  current  (Ig).  The  IBCL  is  calculated  using  argon  ion  mass  until  about  102  minutes  into 
the  test  when  according  to  Eq.  [8.10]  criteria  the  ion  mass  becomes  the  average  mass  of  Cqq  +  Cjq 
ions  <mi>  (Table  8.1  contains  <mj>  definition).  Although  this  selection  is  supported  by  the 
analysis  in  Section  8.2  the  transition  from  Ar  beam  to  mostly  fullerene  beam  is  gradual  and  the 
value  of  IbCL  in  Fig-  8.3c  represents  the  theoretical  limits.  It  is  therefore  seen  that  in  the  second 
half  of  the  test  when  the  fullerene  mass  flow  is  sufficiently  high,  the  measured  ion  beam  at  25  to  30 
mA  is  reasonably  close  to  the  theoretical  limit  of  about  42  mA.  It  should  also  be  noted  that  30  mA 
ion  beam  current  exceeds  that  reported  by  Nakayama(n)  by  nearly  a  factor  of  4. 

Figure  8.3d  shows  the  ion  beam  cost  (eg)  whose  calculations  also  involves  the  two 
different  values  of  ion  mass  (for  argon  and  fullerenes)  discussed  in  the  preceding  paragraph.  As 
seen  the  value  of  eg  in  the  second  half  of  the  test  is  unacceptably  high  ranging  between  1700  to 
2300  volts  with  relatively  stable  value  of  2250  volts  at  the  highest  fullerene  flow  rate. 

Figure  8.3e  shows  the  electrical  and  the  propellant  utilization  efficiencies.  Even  operating 
on  pure  Ar  at  the  beginning  of  the  test  resulted  in  a  low  (~10%)  utilization  efficiency  indicating 
high  losses  in  the  discharge  chamber.  The  stable  values  at  the  highest  fullerene  flow  rate  toward 
the  end  of  the  test  are  70%  and  15%  for  the  electrical  and  utilization  efficiencies,  respectively. 
While  both  values  are  not  acceptable  for  space  propulsion,  the  utilization  efficiency  exceeds  that 

reported  by  Nakayama^11-’  by  a  factor  of  2. 

Figure  8.3f  shows  Isp  and  thrust.  As  seen  in  Table  8.1,  Eq.  20,  the  thrust  is  based  on 
measured  beam  current  and  the  ion  mass  which  was  selected  according  to  the  previously  discussed 
criteria.  The  change  in  ion  mass  from  Ar  to  fullerene  ion  is  responsible  for  the  observed  step 
change  in  the  thrust  level  from  about  1  mN  to  3.5  mN.  These  two  values  should  be  viewed  as 
thrust  limits  while  operating  with  Ar  only  and  mostly  fullerenes  toward  the  end  of  the  test.  Both 
thrust  and  Isp  are  low  even  when  operating  on  pure  Ar  at  the  beginning  of  the  test  indicating  high 
losses  in  the  discharge  chamber  along  with  small  open  fraction  of  our  low  cost  grids. 

Variable  Discharge  Current  and  Intergrid  Voltage  Test  (#071696) 

Figure  8.4  and  8.5  present  the  measured  and  calculated  parameters,  respectively  for  a 
representative  test  where  the  discharge  current  was  varied  from  1.5  to  0.7  Amps  (Fig.  8.4a)  and 
the  beam  voltage  was  varied  between  1000  to  500  volts  while  the  accelerator  voltage  was  held  at 
-500  volts  (Fig.  8.4c).  The  data  are  presented  in  the  same  manner  as  for  the  previous  test  (Fig.  8.2 
and  8.3)  and  all  calculations  were  performed  in  an  identical  manner. 

The  fullerene  mass  flow  presented  in  Fig.  8.5a  was  similar  to  the  previous  test  (see  Fig. 
8.3a)  and  integrated  into  cumulative  value  of  about  8.3  grams.  Pre  and  post  test  weighing  of 
fullerenes  in  the  vaporizer  indicated  that  8.4  grams  were  vaporized  giving  good  closure  and 
reasonable  confidence  in  the  fullerene  flow  rate  calculations.  The  beam  current  shown  in  Fig.  8.4e 
(in  addition  to  the  accelerator  grid  current)  is  a  factor  of  2  to  4  lower  than  the  beam  current  in 
Fig.  8.2e.  This  is  true  even  at  the  beginning  of  the  test  when  operating  on  the  cathode  argon  flow 
only,  yet  the  discharge  current  and  the  grid  voltages  were  the  same  during  that  interval  (43  to 
88  minutes  into  the  test  in  Fig.  8.4a  and  8.2a) 

The  reason  for  this  difference  is  the  significantly  lower  Ar  flow  (0.075  mg/sec  versus  0.15 
mg/sec  in  Fig.  8.2a)  which  also  resulted  in  higher  discharge  voltage  (54  volts  in  the  present  test 
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Fig.  8.4  Discharge  Current,  Voltage,  Argon  Mass  Flow,  Grid  Potentials,  Beam  Current, 
Accelerator  Grid  Current  and  the  Vaporizer  Temperature  were  Measured  and 
Plotted  as  a  Function  of  Time  into  Test  (in  minutes)  for  Test  ID  #071696 
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Fig.  8.5  Fullerene  Mass  Flow  (mdotf).  Cumulative  Fullerene  Mass  Vaporized  (rncummf), 
Ratio  of  Fullerene  Mass  Flow  to  Argon  Mass  Flow  (mdotf/mdotAr),  Ratio  of 
Fullerene  Ion  Mass  Flow  to  Fullerene  Mass  Flow  (mdoti/mdotf),  the  as 
Measured  Ion  Beam  Current  (Ig),  the  Theoretical  Space  Charge  Limited  Beam 
Current  (IbclX  the  Beam  Ion  Cost  (eb),  the  Electrical  (etae)  and  Propellant 
Utilization  Efficiency  (etau),  and  Finally  the  ISp  and  Thrust  are  Plotted  as  a 
Function  of  Test  Time  in  Minutes.  The  above  parameters  were  calculated  from 
the  measured  data  in  Fig.  8.4  using  Equations  in  Table  8.1. 
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relative  to  35  volts  in  Fig.  8.2a).  At  high  fullerene  flow  rates  toward  the  end  of  the  test  the 
discharge  voltages  were  about  the  same  for  both  tests  even  though  the  discharge  currents  were  a 
factor  of  3  different  indicating  that  fullerenes  dominated  the  discharge  conditions.  The  beam 
current  decreased  with  increasing  fullerene  flow  rate,  which  is  an  expected  behavior  given  the 
beam  current  dependence  on  the  ion  mass.  The  beam  current  (Fig.  8.4e)  shows  appropriate 
dependence  to  decreasing  intergrid  voltage  but  the  accelerator  impingement  current  stays  relatively 
constant.  The  beam  ion  cost  (Fig.  8.5d)  the  propellant  utilization  efficiency  (Fig.  8.5e)  and  the 
thrust  are  a  factor  of  2  to  3  worse  than  in  the  previously  discussed  test  indicating  that  the  lower 
discharge  current  (0.7  Amps  versus  1.5  Amps  in  previous  test)  produced  insufficient  ion  density  in 
the  discharge  chamber. 

Test  with  Filament  Cathode  (#073196) 

To  explore  the  role  of  Ar  in  the  discharge  chamber  and  to  compare  our  data  with 
others!7-1 1)  a  filament  cathode  described  in  Section  4.3.1  was  used  in  place  of  the  argon  operated 
hollow  cathode.  The  resulting  data  are  shown  in  Figs.  8.6  and  8.7. 

As  seen  in  Fig.  8.6a  the  discharge  was  started  (after  heating  the  filament  cathode)  on  the 
background  pressure  in  the  tank  which  was  approximately  9  x  10~6  torr.  Because  of  the  low 
pressure  the  discharge  voltage  reached  nearly  180  volts.  With  the  grid  voltages  already  established 
the  beam  current  was  about  1.5  mA.  When  Ar  was  introduced  into  the  discharge  chamber  (at 
about  the  same  rate  as  in  the  previous  test  shown  in  Fig.  8.4b)  the  beam  current  rose  to  about  5.5 
mA  and  the  discharge  voltage  dropped  to  approximately  20  volts.  Argon  flow  was  turned  off  at 
approximately  130  minutes  into  the  test  and  the  beam  current  returned  to  nearly  identical  value  as  at 
the  start  of  the  test.  Argon  flow  was  then  restarted  while  the  fullerene  vaporizer  temperature  (TVap 
in  Fig.  8.6f)  was  rising  and  starting  to  deliver  appreciable  fullerene  flow  into  the  discharge 
chamber  (time  into  test  >  170  min  in  Fig.  8.7a).  At  this  point  the  beam  current  was  decreasing 
essentially  duplicating  the  beam  current  level  and  pattern  observed  with  the  hollow  cathode  in  the 
previous  test  (see  Fig.  8.4e  at  143  min  into  the  test).  The  significant  difference  however,  is  that 
the  discharge  current  in  the  present  filament  cathode  test  was  1.5  Amps  (Fig.  8.6a,  145  to  190 
minutes  interval)  while  in  the  previous  test  the  same  beam  current  was  delivered  at  only  0.5  Amp 
discharge  current.  This  indicates  that  even  with  Ar  in  the  discharge  chamber,  the  filament  cathode 
is  negatively  impacting  the  fullerene  ion  production.  This  is  further  reinforced  by  comparing  the 
beam  currents  delivered  at  the  same  discharge  current  of  0.5  Amps.  In  the  present  filament  cathode 
test  (Fig.  8.6a)  this  occurred  at  about  227  min  into  the  test  with  a  beam  current  of  only  about 
2.5  mA  (Fig.  8.6e)  while  at  nearly  the  same  conditions  in  the  hollow  cathode  test  the  beam  current 
was  about  8  mA  (Fig.  8.4e  at  178  min).  It  required  2.5  Amps  discharge  current  to  produce  8  mA 
beam  current  with  the  filament  cathode  (Fig.  8.6a  and  8.6e  at  227  min). 

The  indisputable  evidence  that  filament  cathode  discharge  is  negatively  impacting  the 
fuilerene  ion  production  came  in  post  test  observation  of  fullerene  remnants  in  the  discharge 
chamber.  While  there  was  always  carbon  film  left  on  the  walls  of  the  chamber  with  hollow 
cathode  tests  (see  Section  7.3),  the  amount  of  carbon  residue  after  the  filament  cathode  tests 
approached  the  amount  evaporated  which  provides  explanation  for  the  low  beam  current. 

8.4  Carbon  Film  Deposition  Results 

Several  samples  of  residue  from  the  discharge  chamber  and  samples  from  films  deposited 
on  substrates  using  the  energetic  fullerene  beam  were  sent  to  SRI  International  for  analysis.  The 
results  are  briefly  discussed  next. 

The  samples  from  the  discharge  chamber  were  collected  from  the  discharge  chamber  wails 
following  several  tests  with  the  hollow  cathode.  As  described  in  Section  7.3  they  were  thin  carbon 
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Fie.  8.6  Discharge  Current,  Voltage,  Argon  Mass  Flow,  Grid  Potentials,  Beam  Current, 
Accelerator  Grid  Current  and  the  Vaporizer  Temperature  were  Measured  and 
Plotted  as  a  Function  of  Time  into  Test  (in  minutes)  for  Test  ID  #073196 
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Fig.  8.7  Fullerene  Mass  Flow  (mdotf).  Cumulative  Fullerene  Mass  Vaporized  (mcummf). 
Ratio  of  Fullerene  Mass  Flow  to  Argon  Mass  Flow  (mdotf/mdotAr),  Ratio  of 
Fullerene  Ion  Mass  Flow  to  Fullerene  Mass  Flow  (mdoti/mdotf),  the  as 
Measured  Ion  Beam  Current  (Ib),  the  Theoretical  Space  Charge  Limited  Beam 
Current  (Ibcl)>  the  Beam  Ion  Cost  (£(,),  the  Electrical  (etae)  and  Propellant 
Utilization  Efficiency  (etau),  and  Finally  the  ISp  and  Thrust  are  Plotted  as  a 
Function  of  Test  Time  in  Minutes.  The  above  parameters  were  calculated  from 
the  measured  data  in  Fig.  8.6  using  Equations  in  Table  8.1. 
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flakes  with  typical  dimensions  of  3  to  4  mm  in  length,  1  to  2  mm  in  width  and  0.1  mm  or  less  in 
thickness.  These  samples  were  examined  using  a  variety  of  techniques  including  RAMAN,  X-ray, 
SALI-MS  and  TEM.  The  mass  spectrometry  show  only  a  very  weak  Cf,o  signal  and  required  a 
strong  desorption  laser  power  to  see  it.  We  can  conclude  from  this  that  the  C^o  has  been  mostly 
destroyed. 

The  Raman  measurements  were  done  by  Wolfgang  Bacsa  at  I  PE  in  Lausanne.  Switzerland 
who  has  done  a  lot  of  Raman  work  on  C^o  samples.  The  spectrum  Bacsa  obtained  on  the  sample 
in  the  range  of  500  to  2000  cm-1  is  shown  in  Fig.  8.8.  There  is  no  evidence  for  diamond  like  or 
sp3  bonding  however,  the  bands  are  broadened  by  disorder.  Bacsa  has  not  identified  all  these 
bands  but  he  concludes  that  the  peaks  in  the  1740  cm'1  are  C-0  bonds  and  the  bands  in  the  1 550- 
1620  region  are  graphene  layers  or  fragments.  The  major  C60  lines  are  either  absent  or  very  weak 
indicating  that  the  Cgo  is  indeed  highly  damaged.  However,  he  concludes  that  the  sample  is  not 
amorphous  because  it  contains  so  many  discrete  bands  which  must  be  due  to  an  abundance  of 
molecular  bound  structures.  Purely  amorphous  carbon  would  give  much  broader  Raman  spectra 
than  that  seen  in  this  sample. 

The  X-ray  diffraction  (Fig.  8.9)  also  shows  that  the  C60  structure  is  not  crystalline  C60, 
however  it  does  show  some  scattering  intensity  in  the  region  of  the  220  and  311  directions.  This 
indicates  that  some  ordering  of  the  carbon  atoms  remains  in  the  sample  but  it  is  nearly  completely 
amorphous. 

For  TEM  analysis  the  sample  was  dispersed  by  sonication  in  alcohol  and  mounted  on  fine 
mesh  copper  TEM  grids.  Only  a  few  photos  were  taken  but  they  clearly  show  some  interesting 
structure  including  a  few  contorted  fullerene  like  tubes  that  have  grown  in  this  sample.  There  also 
seem  to  be  other  particle-like  structures  that  are  not  very  well  resolved.  These  might  be  clusters  of 
damaged  buckyballs  but  they  might  also  be  graphitic  particles.  There  is  insufficient  resolution  to 
see  individual  graphitic  layers  with  this  particular  TEM.  The  growth  of  tubes  is  very  surprising 
under  these  conditions  but  it  is  apparently  a  minor  component  of  the  material. 

Temperature  dependent  resistivity  measurements  done  on  this  sample  are  shown  in 
Fig.  8.10.  These  measurements  appear  appropriate  for  disordered  graphite  however  they  are 
consistent  with  the  results  of  Maiken  and  Taborek(26)  for  hard  amorphous  carbon  film  deposited 
using  fullerene  beam  on  a  heated  substrate. 

The  analysis  of  carbon  films  deposited  on  stainless  steel  and  plastic  (Kapton  tape) 
substrates  using  energetic  fullerene  ion  beams  (-1000  eV)  are  not  complete.  The  only  diagnostic 
tests  completed  are  the  mass  spectrum  analysis  which  detected  no  fullerenes.  Visually,  the  films 
appear  consistent  with  diamond  like  carbon  (DLC). 

8.5  Discussion  &  Recommendations 

The  combination  of  hollow  cathode  (avoiding  hot  surfaces  such  as  the  filament  cathode) 
and  noble  gas  environment  may  be  the  most  likely  reason  for  the  significantly  better  results  with 
the  hollow  cathode  than  with  the  filament  cathode.  Filament  cathode  is  likely  to  fragment 
fullerenes  in  two  different  ways  —  first  is  the  obvious  thermal  fragmentation  which  occurs  at 
800°C  and  above  and  the  second  is  negative  fullerene  creation  by  contact  ionization  where  the 
electronegative  fullerene  removes  an  electron  from  a  low  work  function  surface.  The  negative  ion 
is  then  trapped  in  the  discharge  chamber  until  it  hits  the  walls,  fragments  and  remains  as  a  residue. 
Additionally  it  is  possible  that  the  presence  of  Ar  in  the  discharge  chamber  provides  energy 
distribution/equilibration  path  from  highly  excited  fullerenes  to  neutrals  thus  reducing  fullerene 
fragmentation. 
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As  was  discussed  in  Section  1.2  the  high  carbon  residue  in  the  discharge  chamber  while 
operating  with  filament  cathode  was  also  reported  by  Anderson,  et  al(7)  but  not  by  Nakayama.l11) 
Our  own  experience  from  Phase  I  of  the  program,  where  we  used  filament  cathode  in  quartz 
discharge  chamber,  is  somewhat  inconclusive  because  fragmented  fullerenes  were  found  only  on 
the  filament  cathode  and  on  the  chamber  walls  nearest  to  the  filament  but  most  of  the  chamber  walls 
were  residue  free. 

Overall  the  performance  of  the  thruster  is  disappointing  even  when  operating  with  the 
hollow  cathode.  While  the  initially  predicted  thrust  (see  case  2  in  Table  4.1)  and  beam  current  are 
within  about  a  factor  of  2  of  the  best  measured  values,  the  experimental  beam  ion  cost  is  facto f  of  5 
higher  than  predicted.  Most  importantly,  despite  the  evidence  that  the  beam  was  primarily 
fullerenes  (based  on  presented  data  and  based  on  the  fact  that  deposited  films  were  not  fullerenes 
but  hard  carbon  films  which  requires  at  least  100  eV  impact)  we  have  not  conclusively 
demonstrated  that  the  discharge  and  the  beam  operated  on  fullerenes  in  the  presence  of  Ar.  To 
resolve  these  issues  and  to  bring  the  thruster  to  predicted  performance  we  recommend  the 
following  future  work: 

1)  Use  mass  spectrometer  or  ( ExB  probe)  to  determine  the  beam  composition  (C6o+- 
C60+  +  i  C-/q+,  Ar+,  etc.).  Such  instrumentation  is  extremely  expensive  and  not 
available  to  this  limited  budget  program. 

2)  Upgrade  the  grid  system  to  one  that  is  dished  to  minimize  screen  and  acceleration  grid 
misalignment  at  operating  temperature  and  increase  the  open  area  fraction.  Grids  with 
1470  number  of  holes  (see  Table  4.2)  as  opposed  to  our  low  cost  grids  with  829  holes 
would  significantly  improve  the  thruster  performance.  The  best  alternative  may  be  high 
open  area  carbon/carbon  grids. 

3)  Improve  electron  and  ion  containment  in  the  discharge  chamber  through  shorter 
discharge  chamber  and  better  magnetic  field  utilization. 

4)  Perform  detail  study  on  the  role  of  argon  mixed  with  fullerene  vapor  in  the  discharge 
chamber.  Other  noble  gases  may  be  worthy  of  investigation.  Neon,  because  of  its 
smaller  ionization  cross-section  (which  should,  according  to  Eq.  [8.8]  ensure  that  only 
fullerene  ions  are  produced)  and  He  because  of  its  higher  thermal  conductivity. 
However,  operating  hollow  cathode  on  He  may  be  problematic 

5)  Develop  a  method  for  direct  measurement  of  fullerene  vapor  pressure  and  hence  mass 
flow.  Capacitive  pressure  gauge  made  of  quartz  is  a  possibility.  We  carried  out  its 
initial  development  under  a  BMDO  sponsored  program. (-7)  The  alternative  way  to 
measure  mass  flow  is  quartz  microbalance  which  is  however,  very  difficult  to  use  in  the 
hot  environment  adjacent  to  the  vaporizer. 

6)  Develop  high  temperature  low  pressure  valve  to  turn  the  fullerene  flow  on  and  off  at 
will.  This  would  greatly  simplify  testing. 

7)  Carry  out  experiments  using  pure  C^o  propellant  (as  opposed  to  C60/C70  mix)  to 
determine  which  species  is  more  likely  to  fragment.  Because  C70  is  not  spherical  it  is 
reasonable  to  expect  that  it  is  more  susceptible  to  fragmentation. 
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9.0  CONCLUSIONS 


A  novel  13  cm  fullerene  fueled,  dc  discharge  ion  thruster  system  consisting  of  propellant 
storage  and  vaporizer  integrated  with  the  thruster  was  designed,  constructed  and  tested. 

A  fullerene  vaporization  study  was  carried  out  concurrently  in  a  separate  apparatus  which 
resulted  into  initial  propellant  preparation  specifications.  The  recommended  steps  include  initial 
propellant  drying  at  600  K  in  vacuum  for  10  hours,  followed  by  sublimation  and  condensation  of 
the  fullerenes.  Only  the  condensed  material  should  then  be  used  as  the  propellant.  The  propellant 
should  be  stored  in  dark  containers  under  inert  gas  to  avoid  surface  polymerization  which 
precludes  predictable  vaporization. 

Several  types  of  tests  were  carried  out  with  the  thruster  system  using  (1)  hollow  cathode  in 
the  discharge  chamber  that  was  operated  on  argon  and  (2)  2%  thoriated  tungsten  filament  cathode. 
The  carbon  residue  inside  the  discharge  chamber  while  operating  with  the  argon  fed  hollow 
cathode  was  finite  but  orders  of  magnitude  below  that  observed  when  operating  with  the  filament 
cathode  which  resulted  in  discharge  chamber  carbon  residue  approaching  the  amount  vaporized. 
This  is  attributed  to  (1)  thermal  fragmentation  of  fullerenes  by  the  filament  cathode,  (2)  possible 
negative  fullerene  ion  formation  and  (3)  due  to  unknown  beneficial  effects  of  argon  when  using 
hollow  cathode. 

Using  the  hollow  cathode,  beam  currents  of  the  order  of  30  mA  were  reached  at  total 
intergrid  voltage  of  1.5  kV  yielding  thrust  level  in  excess  of  3  mN.  Low  utilization  efficiency 
ranging  between  8  to  20%  yielded  low  ISp  of  the  order  of  a  few  hundred  seconds.  This 
performance  was  about  a  factor  of  2  below  that  initially  predicted  but  exceeded  by  about  a  factor  of 
2  the  best  data  published  by  others.  Improvements  in  discharge  chamber  design  and  the  grid 
system  are  certain  to  elevate  this  performance.  Exploratory  experiments  related  to  amorphous 
carbon  film  deposition  on  various  substrates  yielded  positive  results  opening  a  possible  terrestrial 
application  for  this  thruster. 

Continuation  of  this  study  toward  promising  and  practical  fullerene  ion  propulsion  with 
possible  dual  (terrestrial)  use  is  strongly  recommended. 
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I .  Introduction 


The  objective  of  this  project  was  to  model  the  ion  thruster  operated  on  a  C-60  and  A r 
propellant  mixture  by  using  ERC’s  particle  simulation  plasma  code.  The  panicle 
simulation  plasma  code,  which  is  based  on  the  panicie-in-cell  (PIC)  method  coupled  with 
the  Monte  Carlo  technique  for  charge-exchange  collisions,  was  developed  by  Peng  of 
ERC  and  his  coworkers  (Ref.  1-3).  It  has  been  used  for  both  two-  and  three-grid  ion 
thruster  simulations.  A  present,  this  simulation  code  can  only  treat  single  species.  As  a 
approximation,  a  two-dimension,  two-grid  and  single  species  model  was  used  to  calculate 
the  electric  potential  field  and  beam  current  with  C-60  propellant  due  to  its  larger  mass 
and  higher  mass  flow  ratio  as  compared  to  Ar  (9:1  ratio).  The  impingement  current  and 
grid  erosion  of  the  accelerator  grid  were  calculated  using  Ar  since  the  cross  section  of  the 
charge  exchange  collision  between  C-60  and  C-60*  is  not  known.  In  order  to  find  the 
equivalent  C-60  neutral  impingement  current  to  the  accelerator  grid  and  the  accurate 
neutral  density  of  Ar  on  each  computational  mesh  point,  simulations  of  the  neutral  flow 
from  the  discharge  chamber  were  conducted  by  using  our  Direct  Simulation  Monte 
Carlo(DSMC)  code  for  both  C-60  and  Ar  gases,  respectively.  Detail  simulation  results  are 
presented  herein. 


II.  Simulation  Results 

1.  Geometry 

Simulations  were  performed  with  two-grid  ion  optics  operated  on  C-60  and  Argon 
propellant  mixture.  The  grid  pattern  was  chosen  to  be  rectangular,  and  the  diameter  of 
the  thruster  is  130  mm.  Three  configurations  were  chosen  for  the  ion  optics  simulations. 
For  all  three  configurations,  the  accelerator  grid  diameter  was  2  mm,  the  center-to-center 
distance  between  adjacent  accelerator  grid  apertures  was  3  mm,  and  the  thickness  of  both 
the  accelerator  and  screen  grids  was  0.5  mm.  Two  different  screen  grid  diameters  were 
used;  2  mm  and  2.5  mm,  respectively.  Two  different  grid  spacings  were  used;  1  mm  and 


0.75  mm,  respectively.  Different  screen  grid  diameter  and  grid  spacing  results  in 
different  electric  potential  field  around  the  grid  region;  therefore,  it  provides  different  ion 
optics  systems.  To  optimize  the  ion  optics  system  for  a  specific  mission,  more 
simulations  are  required. 


All  the  geometry  and  input  parameters  are  listed  in  the  Table  1  and  2. 


It  is  assumed  that  a  neutral  plasma  exists  at  the  upstream  boundary  of  the  computational 
domain,  with  an  ion  temperature  1000  K  and  an  electron  temperature  of  1.5  eV.  The 
background  pressure  of  the  vacuum  chamber  is  assumed  at  8x1 0~5  Torr  and  the  propellant 
is  mixed  as  ninety  percent  C-60  and  ten  percent  Argon. 


Table  1.  Parameters 


diameter  of  thruster 

0.13  m 

discharge  voltage 

30  V 

ion  temperature 

1000  K 

electron  temperature 

1.5  eV 

tank  pressure 

8  x  10'5  Torr 

propellant 

90%  Q0,  10%  Ar 

Table  2.  Grid  configurations 


1 

case  #  1  (mm) 

case  #2  (mm) 

case  #3  (mm) 

diameter  of  accel.  grid 

2.0 

2.0 

2.0 

diameter  of  screen  grid 

2.0 

2.5 

2.5 

grid  hole  pitch 

3.0 

3.0 

3.0 

grid  spacing 

1.0 

1.0 

0.75 

thickness  of  accel.  grid 

0.5 

0.5 

0.5 

thickness  of  screen  grid 

0.5 

0.5 

0.5 

2.  Potential  field  simulation 

Only  C-60  was  used  in  the  ion  optics  simulation.  The  contribution  of  the  Ar  to  the 
potential  field  was  neglected.  This  assumption  introduced  some  small  errors  in  the 
electric  field  simulation  because  Ar  is  much  lighter  than  C-60.  Since  the  cross  section  of 
the  charge  exchange  collision  between  C-60  and  C-60+  is  not  available  and  its  effect  on 
the  potential  field  is  assumed  to  be  small,  the  contribution  of  the  slow  ion  of  C-60  to  the 
potential  field  was  neglected. 

The  ion  optics  were  simulated  by  using  a  two-dimensional  axisymmetric  panicle 
simulation  code  based  on  the  particle-in-cell  (PIC)  method  coupled  with  the  Monte  Carlo 
technique  for  charge-exchange  collisions.  In  the  simulations,  ion  particles  were 
introduced  into  computation  domain  in  the  upstream  region  of  the  screen  grid  and  then 
were  extracted  and  accelerated  through  grid  apertures  to  the  downstream  of  the  thruster 
by  the  electric  potential  field.  The  potential  field  was  calculated  by  solving  Poisson’s 
equation  at  each  computational  mesh  point.  The  calculated  potential  fields  for  three 
different  grid  configurations  are  shown  in  Figs.  1-3  for  intergrid  potentials  of  1500  V, 
2000V,  and  3000V,  respectively.  For  all  three  cases,  the  accelerator  grid  potential  was 
chosen  at  -500  V. 

1500  V  intergrid  potential  cases 

Figs,  la,  b,  and  c  show  the  calculated  electric  potential  fields  under  three  different 
geometry  configurations.  For  comparisons,  the  same  upstream  plasma  number  density 
was  used  in  these  simulations.  The  calculated  beam  currents  for  these  three  grid 
configurations  are  shown  in  Table  3. 
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Table  3.  Comparisons  of  beam  current  for  different  configuration 
at  intergrid  voltage  1500V  (Va=  -500V,  Vs=1000V) 


case  #1 

case  #2 

case  #3 

theoretical  value  (J_max) 

173  mA 

231  mA 

282  mA 

simulated  value 

106  mA 

178  mA 

196  mA 

upstream  plasma  density 

MEllfiUlBSM 

In  Table  3.  the  maximum  theoretical  current  J_max  is  calculated  by, 
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where  iG  is  the  grid  spacing. 

The  simulated  current  is  lower  than  J_max  because  the  simulation  code  experienced  a 
difficulty  in  simulating  the  beam  current  close  to  the  theoretical  maximum  limit.  The 
calculated  perveance  is  about  77%  of  the  theoretical  maximum  value  for  case  2.  In  order 
to  simulate  higher  beam  current,  a  relative  high  plasma  density  (4xl017  m’3  )  was  input 
into  the  upstream  boundary.  This  high  plasma  density  resulted  in  a  concave  potential 
along  the  upstream  equal-potential  contour  of  the  screen  grid  with  a  cusp  toward  the  grid 
region.  The  concave  potential  with  the  cusp  will  prevent  extraction  of  the  high  beam 
current  and  will  affect  the  half  divergence  angle  of  the  ion  beam.  (In  practice,  the 
accelerator  system  may  reach  the  maximum  perveance  condition  without  the  type  of 
difficulty  the  numerical  simulation  code  experienced.) 


In  Table  3,  the  geometry  1  has  the  least  perveance  while  geometry  3  has  the  most.  This  is 
indicated  by  the  electric  potential  field  around  the  upstream  plasma  sheath  region.  In 
Fig. la  (case  1),  the  potential  concave  is  much  deeper  toward  the  screen  grid  because  the 
upstream  plasma  density  (  4xl017  m'3  )  is  relatively  high  for  case  1.  Hence,  the  plasma 
boundary  is  pushed  further  toward  the  grids  for  case  1  than  for  cases  2  (Fig.  lb)  and  3 
(Fig.  lc),  and  the  beam  current  extracted  by  case  1  is  lower  than  case  2  and  case  3.  The 
ion  beams  are  closer  to  the  accelerator  grid  for  configurations  2  and  3  than  for 
configuration  1. 

2000  V  intersrid  potential  cases 

Figs.  2a,  b,  c  shows  the  equal  electric  potential  contours  and  ion  optics  with  a  total 
voltage  of  2000V.  Again,  an  upstream  plasma  density,  6x10 17  m'3,  was  assumed. 
Because  of  the  higher  intergrid  potential,  higher  beam  currents  can  be  achieved  as 
compared  to  the  previous  ones.  Among  the  three  geometric  configurations,  similar  trend 
is  shown  as  in  the  1500  V  cases,  that  is  configuration  1  extracted  less  beam  current  than 
did  configuration  2  and  3.  The  concave  potential  also  exists  along  the  equal  potential 
contour  upstream  of  the  screen  grid  with  a  cusp  toward  the  grid  region.  The  calculated 
beam  currents  for  the  different  configurations  at  the  intergrid  voltage  of  2000  V  are 
shown  in  Table  4. 


Table  4.  Comparisons  of  beam  current  for  different  configuration 
at  intergrid  voltage  2000V  (V a=  -500V,  Vs=1500V) 


case  #1 

case  #2 

case  #3 

theoretical  value  (J_max) 

265  mA 

353  mA 

431  mA 

simulated  value 

147  mA 

265  mA 

296  mA 

upstream  plasma  density 

f iflMliBW 

6.0  x  1017 1/m3 

6 


3000  V  intersrid  potential  cases 


Figs.  3a.  b,  and  c  show  the  equal  electric  potential  contours  and  beam  ions  with  a  total 
voltage  of  3000  V.  Again,  a  plasma  density  of  lxlO18  m'3  was  used  for  these  calculations. 
Higher  total  accelerating  voltage  can  yield  higher  beam  current.  The  calculated  beam 
currents  for  the  different  configurations  at  the  intergrid  voltage  of  3000  V  are  given  in 
Table  5.  The  maximum  beam  current  for  the  3000V  of  the  intergrid  voltage  at  the  same 
configuration  is  much  higher  than  that  for  1500V  and  2000V  of  the  intergrid  voltages. 


Table  5.  Comparisons  of  beam  current  for  different  configuration 
at  intergrid  voltage  3000V  (Va=  -500V,  Vs=2500V) 


case  #1 

case  #2 

case  #3 

theoretical  value  (J_max) 

483  mA 

644  mA 

786  mA 

simulated  value 

274  mA 

471  mA 

502  mA 

UDStream  plasma  density 

1.0  x  1018 1/m3 

1.0  x  1018  1/m3 

From  the  ion  optics  point  of  view,  the  geometry  of  configuration  1  has  poorer 
performance  than  that  of  configurations  2  or  3.  Therefore,  configuration  1  is  not 
recommended  (only  from  the  better  ion  optics  point  of  view).  Between  configurations  2 
and  3,  configuration  3  provides  higher  beam  perveance  while  configuration  2  requires 
less  of  grid  alignment.  Since  the  vacuum  chamber  pressure  is  at  about  8x1  O'3  Torr, 
configuration  3  might  experience  significant  intragrid  arcing  because  of  its  shorter  grid 
spacing.  Therefore,  the  geometry  of  configuration  2  with  total  intergrid  voltage  of  2000 
V  was  selected  for  the  grid  erosion  simulation.  And,  the  calculated  potential  field  of  C- 
60  at  the  intergrid  potential  of  2000V  for  case  2  was  saved  for  the  further  grid  erosion 
simulation  using  Ar. 


3.  Neutral  flow  simulation 


In  order  to  find  the  equivalent  neutral  impingement  current  of  the  accelerator  grid  due  to 
the  atom  C-60  and  the  accurate  neutral  density  of  Ar  on  each  computational  mesh  point, 
the  neutral  flow  lost  from  a  discharge  chamber  was  simulated  by  using  Direct  Simulation 
Monte  Carlo  (DSMC)  type  of  code  for  both  C-60  and  Ar  gases,  respectively.  In  the 
calculations,  the  propellant  was  assumed  to  be  a  pure  gas  and  escaped  from  the  discharge 
chamber. 

The  equivalent  neutral  impingement  currents  of  the  accelerator  grid  (J_na)  were 
calculated  for  C-60  in  two  grid  configurations.  The  calculated  results  are  shown  in  Table 
6.  The  J_na  for  case  2  is  higher  than  that  for  case  1  because  the  screen  grid  hole  diameter 
in  case  2  is  larger  than  the  accel  grid  hole  diameter.  Therefore,  there  are  more  neutral 
particles  striking  the  accel  grid.  Since  the  energy  of  a  neutral  particle  is  much  lower  than 
the  sputtering  threshold,  the  higher  equivalent  neutral  impingement  current  of  the 
accelerator  grid  in  case  2  will  not  cause  much  damage  to  the  grid.  Therefore,  case  2  is 
more  favorable  than  case  1  because  of  its  better  ion  optics. 

Table  6.  Equivalent  neutral  impingement  current  of  the  accel.  grid  due  to  the 
neutral  loss  from  the  discharge  chamber 


case  #1  (mA) 

case  #2  (mA) 

J_na 

21.2 

48.2 

The  neutral  density  distribution  of  the  argon  on  the  computational  mesh  points  was 
calculated  for  the  grid  configuration  2  and  was  used  for  the  calculation  of  charge 
exchange  collisions  and  the  calculation  of  impingement  current  and  erosion  of  the 
accelerator  grid.  Since  the  propellant  utilization  efficiency  is  0.8  and  only  10%  of  the 
propellant  is  the  argon,  the  neutral  density  due  to  neutral  loss  from  the  discharge  chamber 
is  much  less  than  the  neutral  density  caused  by  the  background  pressure  at  8  x  10'5  Torr. 
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Therefore,  the  grid  erosion  of  the  accelerator  grid  is  mainly  determined  by  the 
background  pressure. 

4.  Grid  erosion  simulation 

The  grid  erosion  of  the  accelerator  grid  due  to  slow  charge  exchange  ions  was  calculated 
for  configuration  2  at  an  intergrid  potential  2000V.  Four  charge  exchange  collisions  can 
happen  during  the  extraction  and  acceleration  processes, 

C-60+  (fast)  +  Ar  (slow)  —  C-60  (fast)  +  Ar+  (slow) , 

Ar+  (fast)  +  C-60  (slow)  —  Ar  (fast)  +  C-60+  (slow), 

C-60+  (fast)  +  C-60  (slow)  —  C-60  (fast)  +  C-60+  (slow) , 

Ar+  (fast)  +  Ar  (slow)  - —  Ar  (fast)  +  Ar+  (slow). 

Since  the  first  two  collisions  have  much  less  collision  cross  sections  than  the  last  two  and 
the  charge  exchange  collision  cross  section  for  the  third  one  is  not  available.  Only  the 
fourth  one  was  considered  in  the  simulation.  This  collision  rate  is  determined  by  the 
number  density  of  Ar  and  Ar+  .  Therefore,  the  simulation  of  the  grid  erosion  was 
performed  for  argon  at  10%  of  the  total  flow  rate  under  the  electric  potential  field  from 
the  previous  simulation  for  C-60.  All  parameters  used  for  the  calculations  axe  listed  in 
Table  7.  The  calculated  impingement  current  of  the  accelerator  grid  is  0.3  mA  (Table  3). 
This  is  approximately  1.1%  of  beam  current.  The  calculation  was  under  an  assumption 
that  the  downstream  plasma  density  is  0.025%  of  the  Ar  density  at  the  background 
pressure  8xlO'6Torr.  The  calculated  mass  loss  of  the  molybdenum  accelerator  grid  due  to 
the  erosion  is  8.54x1 0-6  kg  per  100  hours. 


Table  7.  Parameters  for  simulation  of  grid  erosion 


accel.  grid  hole  dia. 

2.0  mm 

screen  grid  hole  dia. 

2.5  mm 

grid  hole  pitch 

3.0  mm 

grid  spacing 

1.0  mm 

thickness  of  accel.  grid 

0.5  mm 

thickness  of  screen  grid 

0.5  mm 

intergrid  voltage 

2000  V 

diameter  of  thruster 

0.13  m 

propellant  utilization 

0.8 

discharge  voltage 

30  V 

ion  temperature 

1000  K 

electron  temperature 

1.5  eV 

tank  pressure 

8  x  10"6  Torr 

propellant 

Ax 

Table  8.  Simulated  results  for  grid  erosion  by  using  parameters  listed  in  the  Table  7. 


beam  current 

26.5  mA 

accel.  impingement  current 

0.3  mA 

half  divergence  angle 

7.5  ° 

erosion 

8.54  mg/lOOhrs 

perveance  per  hole 

1.2  x  1 0‘ 10  A/V3/2 

Hi.  Conclusions 

From  the  ion  optics  point  of  view,  the  geometries  of  configurations  2  and  3  are  more 
favorable  than  that  of  configuration  1.  Because  of  the  high  background  pressure  (8x1  O'5 
Torr),  configuration  3  with  less  intergrid  spacing  might  experience  intragrid  arcing. 
Therefore,  the  configuration  2  was  selected  for  further  grid  erosion  study.  The  reason  to 
chose  an  intergrid  voltage  with  2000  V  is  also  the  result  of  the  high  background  pressure 
concern  or  intragrid  arcing  concern.  The  actual  maximum  total  voltage  without  the 
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intragrid  arcing  under  the  background  pressure  at  8x1 0‘3  can  only  be  found  out  during  the 
practical  thruster  operations. 


Finally,  more  simulations  are  required  to  provide  a  better  definition  of  optimum  ion 
optics  configuration  for  thrusters  operated  with  C-60  and  Ar  mixture  propellant. 
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Configuration:  #1 

Screen  grid  potential:  1000  V 

Accel,  grid  potential:  -500  V 

Upstream  plasma  density:  4.0xl017  m'3 
Beam  current:  106  mA 
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Figure  la.  The  potential  contour  and  beam  ions  for  the  grid  configuration  1  at  the 
intergrid  voltage  1500V. 


Configuration: 

#2 

Screen  grid  potential: 

1000  V 

Accel,  grid  potential: 

-500  V 

Upstream  plasma  density: 

4.0xl017 

Beam  current: 

178  mA 

potential  contours  (Uolt) 
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Figure  lb.  The  potential  contour  and  beam  ions  for  the  grid  configuration  2  at  the 
intergnd  voltage  1500V. 


Configuration:  #3 

Screen  grid  potential:  1000  V 

Accel,  grid  potential:  -500  V 

Upstream  plasma  density:  4.0xl017  m'3 
Beam  current:  196  mA 


Figure  lc.  The  potential  contour  and  beam  ions  for  the  grid  configuration  3  at  the 
intergrid  voltage  1500V. 


Configuration:  #1 

Screen  grid  potential:  1500  V 

Accel,  grid  potential:  -500  V 

Upstream  plasma  density:  6.0xl017  m'3 

Beam  current:  147  mA 
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Figure  2a.  The  potential  contour  and  beam  ions  for  the  grid  configuration  1  at  the 
intergrid  voltage  2000V. 
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Configuration:  #2 

Screen  grid  potential:  1500  V 

AcceL  grid  potential:  -500  V 

Upstream  plasma  density:  6.0xl017  ra‘3 

Beam  current:  260  mA 


Figure  2b.  The  potential  contour  and  beam  ions  for  the  grid  configuration  2  at  the 
intergrid  voltage  2000V. 
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Configuration:  #3 

Screen  grid  potential:  1500  V 

Accel,  grid  potential:  -500  V 

Upstream  plasma  density:  6.0xl017  m'3 
Beam  current:  296  mA 


Figure  2c.  The  potential  contour  and  beam  ions  for  the  grid  configuration  3  at  the 
intergrid  voltage  2000V. 
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Configuration:  #1 

Screen  grid  potential:  2500  V 

Accel,  grid  potential:  -500  V 

Upstream  plasma  density:  l.OxlO18  m‘3 

Beam  current:  274  mA 
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Figure  3a.  The  potential  contour  and  beam  ions  for  the  grid  configuration  1  at  the 
intergrid  voltage  3000V. 
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Configuration: 

#2 

Screen  grid  potential: 

2500  V 

AcceL  grid  potential: 

-500  V 

Upstream  plasma  density: 

l.OxlO18 

Beam  current: 

471mA 

Configuration:  #3 

Screen  grid  potential:  2500  V 

AcceL  grid  potential:  -500  V 

Upstream  plasma  density:  l.OxlO18  m'3 

Beam  current  502  mA 


Figure  3c.  The  potential  contour  and  beam  ions  for  the  grid  configuration  3  at  the 
intergrid  voltage  3000V. 
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APPENDIX  B 


SPREAD  SHEET  ANALYSIS 

SIZING  AND  PERFORMANCE  ESTIMATES  FOR  THE 
FULLERENE  ION  THRUSTER 


The  equations  for  spread  sheet  analysis  of  the  FIT  engine  are  shown  in  Table  B.l.  The 
fullerene  vapor  temperature  in  the  discharge  chamber  (Tn  in  Row  #20)  is  manually  iterated  until  the 
ratio  in  Row  #45  becomes  approximately  1.  During  experiments  the  discharge  chamber 
temperature  must  be  kept  above  Tn  to  avoid  condensation.  The  proper  vaporizer  temperature  (Tv 
in  Row  #59)  is  obtained  by  iteratively  guessing  its  value  until  the  ratio  in  Row  #63  becomes 
approximately  1.  This  means  that  the  required  C60  mass  flow  calculated  in  Row  #52  is  delivered 
by  the  vaporizer  (Row  #62). 
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Table  B1 

Fullerene  Ion  Thrusler/Performance  Equations 


1 

IK1 

CONSTANTS 

Id 

epsO(C2/Nm2)  vac.  permitivity 

0.0000000000088544 

Id 

go(m/sec2)  gravity 

9.81 

7 

qe(C)  elementary  charge 

1.6021E-19 

8 

atomic  unit  mass  amu  (kg) 

!  1 .6604E-27 

k(J/degK)  Boltzman  constant 

i  1 .38E-23 

1 

INPUT/CASE  NO 

T_ 

Dact.(m)  grid  active  dia 

jO.I  3 

ds(m)  screen  grid  hole  dia. 

j0. 002 

da(m)  acc.  grid  hole  dia. 

0.002 

p(m)  grid  hole  pitch 

0.003 

lg(m)  grid  spacing 

10.002 

Vt(Volt)  total  inter,  grid  voltage 

1000 

Vn(volts)  net  voltage=mi*cA2/e 

,500 

eatu(-)  propellant  utilization 

0.8 

Tn  (k)  guess  at  dis.  cham.  Tvapr-ITERATE 

662 

epsbi(volts)  loss  per  beam  ion 

500 

Vdis. (Volts)  guess  discharge  voltage 

30 

etabase(-),  base  eff. 

0.8 

Mw(kg/kgmol)  propellt  mol.  wt. 

720 

dstar  vap(mm)  vaporizer  orifice  dia. 

5.8 

Mwcathode(kg/kgmole)  c.  gas  mole  weight 

39.9 

fcn=(n  of  Ar)/(n  of  C60)  cathode  flow  tract. 

'2.05 

(vary  fen  to  get  at  most  3sccm  cathode  flow) 

OUTPUT -THRUSTER 

fo(-)  open  area  fraction 

=(3.1  4*B1 4A2)/(4*B1  5A2) 

Aat(m2)  surface  area  of  grid 

=3.1 4*B1 2  A2/4 

Aao(m2)  open  area  of  acc.  grid 

r=B30*B31 

N  (-)numbr  of  holes 

=4*B32/(3.14*B14A2) 

mi(kg)  mass  of  prop,  ion 

=B24*$B$8 

ci(m/sec)  ion  exit  speed 

=SQRT  (2*$B$7*B  1 8/B34) 

le(m)  effective  grid  spacing 

=  (B16A2+0.25*B13A2)A0.5 

J(A/m2)  beam  current  density 

=0.62854*SB$5*SQRT(SB$7/B34)*B1 7A1 ,5/B36A2 

Ibeam  (A)  beam  current 

=B37*B32 

ni(f/m3)  ion  number  density 

=B37/($B$7*B35) 

cn(m/sec)  speed  of  neutrals 

=SQRT(2*$B$9*B20/B34) 

ni/nn(-) 

=  (B40/B35)*B1 9/(1  -  B1 9) 

nn(1/m3)  neutral  number  density 

=  B39/B41 

ntot(1/m3)  tot  number  density 

=B42+B39 

ntot(1/m3)  based  on  Tdis.  chamb.  input 
ntot(=R42)/ntot(=R43)iteterate  to  =1 

Tdis.  cavity  (degC)  based  ntot 

=  (5368000000/$BS9)*(1/B20)*EXP(-38000/((8.314/4.186)*B20)) 

=843/B44 

=B20-273 

p(N/m2)  in  dis.  cavity  based  on  ntot 

=B43*SBS9*B20 

p(Torr)  in  dis.  cavity  based  on  ntot 

=0.007451  *B47 

F(N)  thrust  from  eq.4 

=B39*B34*B35A2*B32*(1  -B40/B35+B40/(B35*B1 9)) 

F(grams))  thrust  from  eq.4 

=B49*1000/9.81 

mdot(kg/sec)  mass  flow  eq.5 

=B39*B34*B35*B32/B19 

mdot(mg/sec)  mass  flow  eq.5 

=  B51*1 000000 

Isp(sec)  specific  impls 

=B49/(851  *9.81 ) 

etat(-)  thrust  efficiency 

=B23/(1  +2*$B$7*B21/(B34*B35A2)) 

Pttl=FA2/(2*mdot*eta) 

=  B49A2/(2*B51*B54) 

Pdis=Pttl-Pbeam=Pttl-lbeam*Vnet 

=B55-B38*B1  8 

ldis.(amps)=Pdis/Vdis 

=B56/B22 

OUTPUT -VAPORIZER 

guess  Tv(degC) 

651 

Pc60(mTorr) 

=40000000000*EXP(-38000/((8.31 4/4.1 86)*(B59+273))) 

Pc60(N/m2) 

=  B60/0. 007451/1000 

mdot  c60  (mg/sec) 
mdot/mdot(from  rowf/52) 

[iterate  Tv  to  get  Mdot/mdot=1) 

=  1000000*B61*((3.14*(B25/1 000)  A2)/4)*SQRT  (B  2  4/(2  *3. 1  4*831 4*  (B  59+ 2  73))) 

=B62/B52 

OUTPUT  -  cathode  max  flow 

Nc.gas(1/m3)  cathode  gas  number  densty 

=B27*B43 

cc(m/sec)  thermal  speed  of  cathode  gas 

=SQRT  (2*$B$9*B20/(B26*$B$8)) 

mdot.cgas(mg/sec) 

=1 000000*B66*B26*$BS8*B67*B32 

Vdot  cgas(sccm)  (keep  at  min.=about  2sccm) 

=60*B68/(1 02000*B26/(831 4*293)) 

mass  fration  (mdot  cqas/mdotC60) 

=B68/B52 
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